FILE  CORy 

U  ☆ 


ACTIVE  VIBRATION  CONTROL  OF  A  CANTILEVERED 
BEAM  WITH  THREE  ELASTIC  COORDINATES 
THESIS 

Thomas  A.  Cristler 
Captain,  USAF 

AFIT/GA/AA/87D-1 


CT !  nr 


FEB  0  8  1953 


DEPARTMENT  OF  THE  AIR  FORCE 

AIR  UNIVERSITY 

AIR  FORCE  INSTITUTE  OF  TECHNOLOGY 


Wright-Patterson  Air  Force  Base,  Ohio 


88  2  OS  064 


v-v'>.‘-  r-'.v  .N  r- 


AFIT/GA/AA/87D-1 


ACTIVE  VIBRATION  CONTROL  OF  A  CANTILEVERED 
BEAM  WITH  THREE  ELASTIC  COORDINATES 
THESIS 


Thomas  A.  Cristler 
Captain,  USAF 

AFIT/GA/AA/87D-I 


Approved  for  public  release;  distribution  unlimited 


DT!C 

ELF'"."': 

FEB  0  8  i988 

H 


AFIT/GA/AA/87D-1 


ACTIVE  VIBRATION  CONTROL  OF  A  CANTILEVERED  BEAM  WITH  THREE  EL<^STIC 

COORDINATES 


THESIS 


Presented  to  the  Faculty  of  the  School  of  Engineering 
of  the  Air  Force  Institute  of  Technology 
Air  University 

In  Partial  Fulfillment  of  the 
Requirements  for  the  Degree  of 
Master  of  Science  in  Astronautical  Engineering 


Thomas  A.  Cristler,  B.S. 
Captain,  USAF 


^  Accession  For 

NT  IS  GFA.S-.I 
DTIC  TA!'. 

j  UllH!  ilio  •  .  ■'! 

j  \f.-:  •  '  .  1 


December  1987 


Approved  for  public  release;  distribution  unlimited 


Preface 


This  experiment  was  developed  and  run  in  conjuction  with  the  Flight 
Dynamics  Laboratory,  Vibration  and  Acoustics  Branch,  Wright  Patterson 
AFB.  Ohio,  as  a  part  of  their  in-house  Large  Space  Structures  Technology 
Program,  The  experiment  was  conceived  and  constructed  in  Building  24C. 
Area  B,  Wright  Patterson  AFB.  Project  Manager  for  the  experiment  is  Mr. 
Robert  Gordon  (AV785 - 52 36 ) ,  and  the  Lead  Engineer  is  Mr.  Wayne  Yuen. 

The  purpose  here  is  not  only  to  document  the  work  done  to  develop 
and  test  the  system,  but  to  provide  those  future  researchers  interested 
it.  conducting  experiments  with  this  aparatus  with  a  working  level 
knowledge  of  the  system,  its  development,  and  some  of  its 
i ri i os vnc rac ies .  WTiile  the  text  may  seem  to  get  bogged  down  with 
incredible  detail  on  actuator  development  and  structural  modelling, 
tliere  are  some  issues  raised  about  these  subjects  that  even  the  reader 
interested  only  in  Large  Space  Structure  Control  applications  should 
fir-t'.  useful.  However,  unless  one  is  planning  on  using  the  experiment 
for  future  investigations,  studying  Chapters  2  and  3  in  great  detail  is 
T  •  Tiecessarv  for  understanding  the  subsequent  material. 

F  ,ic i 1 i ar i za t i on  with  the  system  configuration  and  problems  encountered 
ii-  df'.-eloping  it  are  probably  the  most  important  points  to  be  gleaned 
fjorr  those  sections  of  this  report. 

Uliile  a  number  of  people  made  significant  contributions  to  this 
project,  I  would  like  to  thank  several  individuals  in  particular  for 
tFieir  help,  understanding,  and  sympathetic  listening  when  things  weren't 
working  as  planned.  First,  my  advisor.  Dr.  Robert  Calico  for  guiding  me 
through  this  often  frustrating  project.  Bob  Gordon  and  Wayne  Yuen 
provided  me  not  only  with  sound  technical  advice,  but  did  most  of  the 
dirty  work  of  procuring  hardware  and  lining  up  technical  services  and 


support.  And  of  course,  the  people  who  really  make  an  experiment  work, 
no  matter  how  much  credit  the  engineers  may  claim.  Earl  Rodgers  and  Mike 
Banford,  the  lab  technicians.  Finally,  a  special  thanks  to  my  friend 
Uendy  Motlong,  and  the  members  of  GA-87D,  all  of  whom  made  this  last 
year  and  a  half  a  bearable  experience. 


Thomas  A.  Cristler 
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Abstract 


Many  control  methods  have  been  proposed  for  dealing  with  the  large 
space  structure  vibration  control  problem.  To  experimentally  evaluate 
these  various  approaches  in  a  way  which  will  allow  consistent  comparison 
of  results  requires  a  baseline  experiment  in  which  all  variables  are 
understood  and  controlled.  From  this  baseline,  the  various  aspects  of 
each  coTitrol  scheme  can  be  implemented,  and  their  relative  merits 
compared  on  a  consistent  basis.  This  experiment  was  implemented  using  a 
vertically  suspended  cantilever  beam  with  rectangular  cross  section. 
Proof  mass  actuators  were  developed  to  provide  control  force  inputs  to 
the  structure.  Closed  loop  control  was  formulated  using  linear 
quadratic  regulator  theory  and  results  are  compared  with  simulation  and 
eigen^’alue  predictions  to  establish  baseline  performance.  Modal 
suppression  techniques  were  implemented  to  demonstrate  control  of 
selected  modes  while  maintaining  overall  system  stability.  Results 
applicable  to  future  testing  and  development  in  the  large  space 
structure  control  area  are  identified. 
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ACTIVE  VIBRATION  CONTROL  OF  A  CANTILEVERED  BEAM  WITH  THREE  ELASTIC 


COORDINATES 


I .  Introduction 


As  man's  presence  in  space  grows,  so  does  the  size  of  the  vehicles 
he  deems  necessary  to  put  there.  Programs  such  as  the  Space  Defense 
Initiative  and  NASA's  space  station  are  generating  requirements  for  very 
large  vehicles  which  must  point  with  arc-second  accuracies,  hold  the 
shape  of  optics  surfaces  to  within  a  few  tenths  of  the  electromagnetic 
wavelength  being  used,  or  maintain  a  vibration  free  environment  for 
experimentation  while  docking  and  cargo  transfer  occurs  elsewhere  on  the 
vehicle.  While  these  operations  are  readily,  and  in  some  cases  commonly 
accomplished  on  earth,  they  present  a  challenging  problem  when  attempted 
in  the  space  environment. 

As  platforms  grow  in  size,  economics  dictate  a  reduction  in  weight. 
The  result  is  a  large  flexible  structure  which  tends  to  have  many 
closely  spaced,  lightly  damped  vibration  modes  that  are  easily  excited. 
Mechanisms  for  inducing  vibration  in  the  structure  range  from  the 
attitude  control  system  and  active  power  generation,  to  cargo  transfer, 
or  personnel  moving  about  the  vehicle.  If  the  amplitude  or  duration  of 
the  vibration  will  degrade  the  performance  of  the  system,  then  some 
means  of  controlling  it  must  be  introduced.  Adding  passive  damping 
(shock  absorbers,  visco-elastic  materials,  etc.)  would  most  likely  be 
the  first  attempt  at  a  solution.  However,  there  are  cases  where  passive 


damping  will  still  not  achieve  the  desired  performance  and  an  active 
means  of  control  must  be  introduced. 

Many  approaches  have  been  proposed  for  dealing  with  the  active 
vibration  control  problem,  and  there  is  considerable  ongoing  work  in  the 
field.  Two  examples  of  current  work  include  Lockheed's  demonstration  of 
a  high  authority  control/low  authority  control  (HAC/LAC)  scheme  on  a 
suspended  plate,  and  static  and  dynamic  shape  control  of  a  fixed- free 
beam,  which  was  accomplished  at  the  Jet  Propulsion  Laboratory  using 
Kalman  and  adaptive  control  techniques  (1:3).  Many  of  the  experiments 
conducted  to  date  have  used  ground  referenced  sensors  or  actuators, 
which  is  not  representative  of  the  problem.  While  this  simplifies  the 
implementation  of  the  experiment,  it  seems  to  be  oriented  more  toward 
control  law  evaluation  than  being  a  complete  emulation  of  the  large 
space  structure  control  problem. 

To  provide  a  more  realistic  test  bed  for  research  in  this  area.  Air 
Force  Wright  Aeronautical  Laboratory,  Wright-Patterson  AFB,  initiated  an 
in  house  program  under  its  vibrations  branch  (AFWAL/FIBG)  to  investigate 
large  space  structures  technologies.  An  advanced  beam  experiment  was 
devised  to  demonstrate  active  vibration  control  of  a  cantilevered  beam 
in  two  orthogonal  bending  axes  as  well  as  torsion  about  its  long  axis. 
This  configuration  has  closely  spaced  and  coupled  modes  and  features 
inertial  sensors  and  actuators.  The  experiment  was  designed  to  emulated 
as  many  large  space  structure  problems  as  possible.  AFIT  was  invited  to 
participate  in  the  implementation  and  testing  of  this  experiment,  which 
is  depicted  in  Figure  1.1. 

Two  related  methods  of  structural  control  which  have  received 
attention  in  recent  years  are  reduced  order  and  decoupled  controllers. 
The  importance  of  these  controllers  lies  in  the  fact  that  a  large  space 
structure  will  in  general  present  many  more  modes  than  can  reasonably  be 


controlled  by  a  single  system.  Weight  and  hardware  capability  become 


serious  issues  when  the  modes  to  be  controlled  can  number  in  the 


hundreds.  If  the  system  can  be  subdivided  and  a  subset  of  the  modes  to 


be  controlled  assigned  to  separate  controllers,  the  implementation  of 


the  structural  control  can  be  greatly  simplified.  Indeed,  not  all  of 


the  modes  may  require  active  control  to  meet  the  vehicle  performance 


requirements.  In  such  cases,  reducing  the  order  of  the  control  model 


can  greatly  reduce  the  size  and  weight  of  the  requisite  control  system. 


The  difficulty  with  implementing  this  reduced  order  model  is  that 


the  sensors  may  still  contain  information  about  the  modes  which  have 


been  omitted  from  the  model.  Also,  since  the  actuators  operate  on  a 


continuous  structure  they  will  excite,  and  potentially  destabilize  these 


ignored  modes.  These  effects  have  been  classified  as  observation  and 


control  spillover  respectively  (Ref  2).  Work  by  Coradetti  (Ref  3) 


showed  that  spillover  could  be  eliminated  by  finding  a  transformation 


matrix  which  is  applied  to  the  feedback  gains,  thus  "suppressing"  the 


omitted  modes.  Calico  and  Janiszewski  (Ref  4)  showed  that  eliminating 


either  observation  or  control  spillover  was  sufficient  to  ensure 


stability  of  the  suppressed  modes  and  demonstrated  a  procedure  for 


calculating  the  appropriate  transformation  matrix,  Wright  (Ref  5) 


implemented  this  modal  suppression  technique  on  a  single  bending  axis  of 


a  cantilevered  beam.  While  being  able  to  show  increased  stablilitv  in 


the  second  bending  mode  of  the  beam  without  destabilizing  the  first  or 


third  modes,  direct  correlation  between  predicted  and  measured 


performance  was  not  achieved. 


The  goals  of  the  research  presented  here  are  to  conduct  the  system 


identification  of  the  experiment,  to  include  the  actuators,  sensors,  and 


structure,  and  demonstrate  predictable  closed  loop  control  response 


using  state  space  control  techniques.  Once  this  has  been  accomplished. 


^  m  •**.  *■"» 


implementation  and  demonstration  of  reduced  order  controllers  will  be 
attempted,  leading  up  to  a  demonstration  of  a  decoupled  control 
algorithm . 


II.  Ac  tuators 


Actuator  System  Description 

The  actuators  chosen  by  AFWAL  for  use  on  this  system  are  proof 
mass  linear  motors  which  provide  control  force  using  momentum  exchange 
between  the  base  and  the  moving  mass.  The  specific  actuators  being  used 
in  this  experiment  are  based  on  a  TRW  design  used  in  the  AFWAL  sponsored 
VCOSS  program.  System  drawings  and  descriptions  were  provided  to  AFWAL, 
who  purchased  and  fabricated  the  necessary  components  to  assemble  the 
devices . 

The  actuators  consist  of  a  linear  motor  coil  mounted  on  two  support 
brackets  (see  Figure  2.1).  A  cylindrical  proof  mass  of  0.9  kg  contains 
the  motor  magnets  and  is  driven  by  the  motor  coil.  This  mass  travels  on 
linear  bearings  along  a  center  shaft  with  a  nominal  travel  of  ±0.5 
inches.  The  nominal  motor  parameters  published  by  the  manufacturer, 
Kaiser  Electroprecision,  are  listed  in  Appendix  A.  The  actuator  is 
instrumented  with  a  Linear  Variable  Differential  Transformer  (LVDT) , 
manufactured  by  Schaevitz  Engineering,  which  provides  feedback  of 
relative  position  between  the  proof  mass  and  the  motor  base.  The  proof 
mass  itself  is  instrumented  with  an  Endevco  piezoresistive  accelerometer 
to  provide  feedback  control  of  the  proof  mass  acceleration.  LVDT  and 
accelerometer  specifications  and  calibrations  are  listed  in  Appendix  A. 
The  motor  coil  is  driven  by  a  power  amplifier  circuit  (see  Appendix  A) 
which  transforms  a  voltage  command  into  a  drive  current.  The  power 
amplifier  has  a  current  limiter  to  prevent  burning  out  the  motor  coil 
which  has  a  steady  state  current  limit  is  2  amps.  This  limits  the  force 
output  of  the  actuators  to  h  Ibf. 
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Initial  open  loop  testing  of  the  actuators  found  them  to  have 


several  undesirable  characteristics: 


1)  The  low  frequency  output  was  very  non-linear.  The  bearing 


friction  and  hysteresis  in  the  system  caused  the  proof  mass  to  wander 


away  from  the  center  of  travel  when  commanded  at  low  frequencies.  Since 


the  actuator  stroke  is  limited,  the  amplitude  of  low  frequency  commands 


must  be  kept  small  to  prevent  driving  the  proof  mass  into  the  support 


brackets.  However,  when  the  input  is  kept  low,  friction  can  overcome 


the  proof  mass  motion  and  cause  it  to  stick,  or  at  least  signif icantly 


lag  the  sinusoidal  input  command.  This  leaves  the  proof  mass  off 


center,  and  after  a  few  cycles  the  mass  is  driven  into  the  end  stop. 


2)  The  zero  reference  potentiometer  in  the  power  amplifer,  which 


provides  centering  adjustment  for  the  proof  mass,  required  continual 


adjustment  during  the  initial  testing  of  the  actuators.  Changes  in  input 


gain  and  changes  in  frequency  would  cause  a  shift  in  the  zero  reference 


(the  center  point  of  the  actuator  travel).  This  problem  was  evident  at 


higher  frequencies  (greater  than  3  Hz),  and  is  a  separate  condition  from 


the  "wander"  previously  described. 


3)  The  actuator  frequency  response  has  a  low  frequency  roll-off 


and  phase  shift  in  the  vicinity  of  the  fundamental  bending  modes  of  the 


structure . 


To  modify  the  frequency  response  of  the  actuators  and  maintain 


proof  mass  centering,  a  feedback  control  system  using  classical  analog 


control  techniques  was  implemented  around  each  actuator. 


Feedback  Compensation  Design  Goals 


The  actuator  control  was  developed  around  accomplishing  the 


following  objectives: 


%•* 
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1)  Actuator  frequency  response  should  be  "flat”  over  the  structure 
control  bandwidth,  i.e.  constant  magnitude  and  zero  phase.  The  control 
bandwidth  for  this  experiment  has  been  limited  to  0  to  50  Hz. 

2)  The  actuator  proof  mass  should  maintain  an  inertial  position 
when  the  actuator  is  not  being  commanded,  which  prevents  the  actuator 
proof  masses  from  contributing  to  the  structure  inertia. 

3)  The  proof  mass  should  maintain  its  centering  when  commanded  at 
different  frequencies  and  amplitudes. 

The  proof  mass  travel  should  be  limited  to  the  actuator  stroke 
of  ±0.5  inches . 
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Open  Loop  Testing, 

Each  actuator  was  tested  ii.  the  original  configuration  to  determine 
its  frequency  response  characteristics.  These  tests  were  performed  by 
driving  the  actuator  power  amplifier  with  a  random  input  signal  and 
measuring  the  acceleration  of  the  proof  mass  with  the  actuator  base 
fixed  to  a  clamped  plate.  The  frequency  response  was  calculated  using 
Fast  Fourier  Transforms  (FFTs)  on  the  measured  data  and  dividing  the 
output  spectrum  by  the  input  spectrum  to  determine  the  transfer 
function.  The  data  collection,  FFT  and  transfer  function  calculation 
w..ie  all  automatically  performed  by  an  Ono  Sokki  CF-910  Dual  Channel  FFT 
Analvzer.  The  output  signal  measured  was  the  voltage  3ut  of  the  proof 
mass  accelerometer.  This  output,  and  consequently  the  transfer 
function,  can  easily  be  transformed  to  force  by  adding  4.9  dB  to  the 
magnitude  response  function  (accelerometer  scale  factor  time  mass  of  the 
proof  mass) . 
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The  transfer  functions  measured  for  each  actuator  are  shown  in 
Figures  2.2  to  2.5.  The  response  is  consistent  between  the  four 
actuators,  with  the  exception  of  the  phase  response  of  actuator  #4,  The 
180°  phase  shift  was  found  to  be  due  to  the  motor  coil  leads  being  wired 
opposite  of  the  other  three  actuators.  When  corrected,  the  phase 
response  of  actuator  #4  was  consistent  with  actuators  #1  through  #3. 

System  Model 

Initially  the  actuators  were  modelled  using  standard  linear  motor 
equations  of  motion.  This  approach  did  not  prove  fruitful,  primarily 
because  the  power  amplifier  circuit  was  constructed  so  that  it 
significantly  reduced  the  back-emf  of  the  motor.  Because  of  this  the 
basic  model  was  developed  by  fitting  the  measured  transfer  function  with 
the  circuit  dynamics  of  the  power  amplifier  plus  any  additional  dynamics 
necessary  to  match  the  response.  From  the  power  amplifier  circuit 
diagram  in  Appendix  A,  the  power  amplifier  frequency  domain  dynamics  are 
found  to  be 

I(s)  500 

E(s)  S  +  2500 

where  I(s)  is  the  current  signal  to  the  motor  coil,  and  E(s)  is  the 
voltage  signal  to  the  power  amplifier. 

The  resulting  analytic  transfer  function  that  gives  the  best  fit 
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Figure  2.5.  Actuator  #4  Open  Loop  Transfer  Function 
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Figures  2.6  and  2.7  show  comparisons  of  the  magnitude  and  phase  response 
of  this  transfer  function  with  the  frequency  response  measured  for 
actuator  #1.  The  model  matches  the  data  very  closely  except  in  the  very 
low  frequencies  (below  1  Hz) .  The  low  frequency  data  is  somewhat 
suspect  due  to  the  short  sample  times  and  low  input  amplitude  required 
to  complete  a  test  on  the  uncompenstated  actuator  without  driving  the 
proof  mass  into  the  end  bracket.  Figure  2.8  shows  the  coherence  for  the 
actuator  #1  frequency  response,  which  is  a  measure  of  the  correlation 
between  the  input  and  output.  A  coherence  of  1.0  indicates  perfect 
correlation  between  input  and  output,  while  0.0  indicates  no  correlation 
at  all.  The  coherence  measured  for  low  frequencies  is  very  poor, 
indicating  the  data  is  suspect. 

The  model  also  does  not  include  the  dynamics  which  are  obviously 
present  in  the  measured  data  around  150  Hz.  This  effect  appears  to  be  a 
mechanical  response  in  the  actuator,  however  the  frequency  is  much  too 
low  to  match  either  the  first  bending  mode  of  the  end  bracket  or  the 
first  axial  vibration  mode  of  the  motor  shaft.  To  isolate  the  effect 
the  mounting  bracket  was  instrumented  with  an  accelerometer  and  the 
response  measured  while  driving  the  actuator  with  a  random  input  signal. 
Figure  2.10  shows  the  transfer  function  between  the  voltage  command 
input  and  the  bracket  acceleration.  Comparing  this  response  with  the 
actuator  transfer  function  indicates  the  energy  is  indeed  being 
transfered  to  the  mounting  brackets  and  not  the  proof  mass  in  the  150  Hz 
range.  Since  this  effect  is  outside  the  desired  control  bandwidth  for 
the  structure,  the  command  signals  should  not  have  any  spectral  content 
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Figure  2.6.  Actuator  »1  Magnitude  Response  -  Measured  vs  Model 
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in  this  frequency  range.  Ignoring  it  should  not  cause  a  problem  other 
than  perhaps  adding  some  noise  in  the  acceleration  feedback. 


To  further  evaluate  this  model,  and  determine  that  the  high 
frequency  response  is  not  aggravated  by  changing  the  system  dynamics,  a 
position  feedback  loop  was  implemented  using  the  LVDT  output  signal. 

The  LVDT  is  modelled  as  a  double  integrator  of  the  proof  mass 
acceleration  and  a' gain.  The  gain  is  calculated  by  applying  the  nominal 
motor  parameters  and  the  nominal  LVDT  scale  factor  (see  Appendix  A) . 

The  resulting  closed  loop  model  is  represented  by  the  block  diagram  in 
F i gure  2.10. 
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Figure  2.10.  Position  Feedback  Block  Diagram 


Figures  2,11  and  2.12  compare  the  transfer  function  derived  from 
tiii.s  analytic  model  and  the  transfer  function  measured  for  actuator  »1  , 
Agreement  between  the  model  and  measured  data  is  good  and  the  high 
frequency  re.sponse  has  not  been  aggrevated. 
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Figure  2.11.  Position  Feedback  Magnitude  Response  Comparison 
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Figure  2.12.  Position  Feedback  Phase  Response  Comparison 
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Control  Design 


The  parameters  accessible  in  this  system  are  the  current  command 
from  the  power  amplifier  to  the  motor  coil,  the  proof  mass  acceleration, 
and  the  relative  position  between  the  end  bracket  and  proof  mass  as 
measured  by  the  LVDT.  Evaluating  the  affect  of  these  parameters  as 
feedback  signals  on  a  heuristic  level,  one  notes  that  position  feedback 
has  the  same  effect  as  adding  a  spring  between  the  end  bracket  and  proof 
mass.  This  would  help  in  accomplishing  the  design  goals  of  maintaining 
the  centering  of  the  proof  mass,  as  well  as  helping  to  limit  its  travel. 
Since  proof  mass  acceleration  varies  from  the  force  by  only  a  constant, 
it  would  have  the  effect  of  feeding  back  the  desired  output.  This 
should  provide  not  only  a  more  accurate  response,  but  also  an  error 
signal  which  will  accomplish  the  goal  of  maintaining  an  inertial  proof 
mass  position  when  the  actuator  is  uncommanded.  A  current  feedback  loop 
around  the  power  amplifier  may  help  improve  the  accuracy  of  the  current 
command,  but  since  the  power  amplifier  dynamics  are  already  outside  of 
the  control  bandwidth,  it  should  not  be  necessary  to  modify  its 
response . 

To  move  the  actuator  dynamics  outside  the  control  bandwidth,  a  root 
locus  analysis  was  used  to  determine  the  effects  of  closing  the  feedback 
loops  discussed  above.  Closing  the  position  feedback  loop  adds  a  low 
order  pole  to  the  system,  which  with  increasing  gain  combines  with  the 
(S  +  35)  root  and  forms  a  complex  pair.  The  break  frequency  for  this 
pair  occurs  at  approximately  6.5  Hz.  If  the  acceleration  loop  is  closed 
around  the  position  feedback,  the  order  of  the  system  remains  the  same 
and  the  break  frequency  of  the  complex  pole  pair  is  lowered  as  the  gain 
is  increased.  To  have  all  effects  of  these  roots  out  of  the  control 
bandwidth  requires  placing  them  below  0.2  Hz. 


Figure  2.13  shows  a  block  diagram  of  this  preliminary 
conf igurat ion .  While  it  is  theoretically  possible  to  move  the  complex 
pole  pair  as  close  to  the  origin  as  desired,  the  gain  required  to  move 
the  poles  below  0.2  Hz  is  approximately  100.  Since  nominal  gain  for  the 
operational  amplifiers  intended  to  be  used  in  this  control  circuit  is 
10,  this  would  require  more  than  one  amplification  stage,  and  would  most 
likely  saturate  the  amplifiers.  An  alternative  to  this  is  to  boost  the 
acceleration  gaiti  as  much  as  the  circuit  will  handle,  and  then  adjust 
the  overall  response  by  prefiltering  the  input  signal.  Applying  a  gain 
of  lU  to  the  acceleration  feedback  would  move  the  poles  below  1  Hz. 
Increasing  the  position  feedback  gain  will  drive  the  damping  ratio  down 
and  reduce  the  frequency  range  over  which  the  phase  contribution  of 
these  poles  acts.  A  predicted  frequency  response  for  this  configuration 
is  shown  in  Figure  2.14.  This  configuration  has  a  prefilter  (lag 
compensation),  position  feedback  gain  of  0.1,  and  an  acceleration 
fetdback  gain  of  10.  This  provides  not  only  a  fairly  flat  response  over 
the  bandwidth,  but  also  the  centering  and  inertial  position  control 
des i red . 

Cl osed  Loop  Testing 

Many  technical  problems  were  discovered  while  implementing  the 
control  design.  A  discussion  of  each  of  the  major  problems  encountered 
and  their  resolution,  if  any,  is  in  Appendix  B.  There  were  two 
important  problems  which  led  to  changing  the  control  design.  First,  the 
acceleration  gain  could  not  be  increased  beyond  a  factor  of  6  without 
high  frequency  noise  dominating  the  feedback  signal.  Figure  2.15  shows 
the  frequency  response  measured  for  the  baseline  configuration  with  the 
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acceleration  feedback  gain  at  5  and  the  position  feedback  gain  at  0.05. 
While  the  magnitude  response  is  fairly  flat,  there  is  too  much  low 
frequency  phase  shift.  The  input  prefilter  can  be  adjusted  to  further 
reduce  the  phase,  but  it  then  has  the  undesireable  side  effect  of 
increasing  the  magnitude  over  the  same  range.  Since  noise  was  the 
limiting  factor  a  low-pass  filter  with  a  break  frequency  of  100  Hz  was 
placed  in  the  acceleration  feedback  loop. 

The  second  problem  was  that  the  centering  of  the  proof  mass  shifted 
as  the  amplitude  of  the  input  signal  was  increased,  even  with  the 
position  feedback  loop.  In  an  attempt  to  correct  the  centering  problem, 
a  feedback  loop  was  placed  around  the  power  amplifier.  Initially,  this 
loop  was  also  plagued  by  high  frequency  noise  problems.  By  changing  the 
inverter  in  the  feedback  path  to  a  low-pass  filter  with  a  break 
frequency  of  50  Hz,  the  noise  problem  was  significantly  reduced.  This 
modification  did  not  help  the  centering  problem,  however  it  did  reduce 
the  overall  noise  in  the  system  and  improve  the  "quality"  of  the 
response  data  (measured  transfer  functions  were  much  smoother  and  more 
repeatable).  The  solution  to  the  centering  problem  was  later  found  to 
be  associated  with  unstable,  very  high  frequency  oscillations  (kilohertz 
range)  in  the  power  amplifier.  Details  of  this  problem  are  discussed  in 
Appendix  B. 

Making  the  modifications  discussed  above  and  several  adjustments  to 
gains  results  in  the  configuration  shown  in  Figure  2.16.  Figure  2.17  is 
the  circuit  implementation  of  this  block  diagram.  Since  the  low-pass 
filters  in  the  feedback  loops  add  zeros  in  the  closed  loop  system,  an 
additional  pole  was  added  as  a  prefilter  to  reduce  the  phase.  The 
desired  response  for  each  actuator  was  set  by  adjusting  the  various 
potentiometers  in  the  circuit,  which  accounted  for  variations  in  gain 
and  sensitivity  of  the  various  elements  in  the  system. 
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The  predicted  frequency  response  of  this  configuration  is  shovm  in 
Figure  2.18.  The  measured  response  of  the  four  actuators  is  overlaid  on 
the  prediction  in  Figures  2.19  and  2.20,  and  show  excellent  agreement 
with  the  prediction.  The  frequency  response  was  measured  using  a  sine 
sweep  rather  than  the  random  input  tests.  This  test  drives  the  actuator 
at  a  single  frequency,  measures  the  response,  and  then  increments  to  the 
next  frequency.  The  method  gives  more  consistent  results  and  is  a  much 
better  approximation  to  the  actual  commands  the  actuator  will  receive 
when  mounted  on  the  structure. 


Characterization  of  Final  Actuator  Conf iRuration 

To  fully  evaluate  the  actuator  performance,  a  series  of  tests  were 
performed  that  simulated  the  operating  environment  an  actuator  will 
experience  when  mounted  on  the  beam.  These  tests  measured  the  actuator 
frequency  response  with  the  base  free  to  move,  the  maximum  force  output 
of  the  actuators  at  each  modal  frequency,  the  uncommanded  actuator 
response  to  base  motion  (how  well  the  proof  masses  "float"),  and  the 
variation  of  the  transfer  function  while  forces  are  applied  orthogonal 
to  the  direction  of  proof  mass  motion. 

A  better  estimate  of  the  actual  performance  the  actuators  will 
deliver  when  fixed  to  a  moving  beam  can  be  made  by  testing  each  actuator 
on  a  slip  table.  The  slip  table  is  a  plate  which  rides  on  bearings  and 
is  free  to  move  in  a  single  direction.  The  test  configuration  is 
depicted  in  Figure  2.21.  The  force  output  of  the  actuator  was  derived 
by  measuring  the  acceleration  of  the  base  plate  and  multiplying  by  the 
combined  mass  of  the  plate  and  actuator  base  (all  components  except  the 
proof  mass).  Figures  2.22  to  2.25  show  the  transfer  functions  measured 
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Compensated  Actuator  Transfer  Function  -  Predicted  vs  Measured  Phase 


Figure  2.23.  Actuator  #2  Slip  Table  Transfer  Function 


Table  Transfer  Func 


for  each  actuator.  For  most  of  the  tests  there  was  very  little  change^ 
from  the  fixed  base  tests,  mainly  a  phase  shift  of  a  few  degrees. 
However,  for  very  low  amplitude  tests  the  response  sometimes  changed 
noticably,  especially  in  magnitude.  This  seemed  to  be  due  to  the  wires 
attached  to  the  actuator  preventing  free  motion  of  the  slip  table.  The 
motor  power  lead  is  a  relatively  stiff  cable,  and  even  though  it  was 
suspended  vertically  from  the  actuator,  it  often  had  more  effect  on  slip 
table  motion  than  the  actuator  itself. 

The  maximum  force  output  of  a  proof  mass  actuator  is  limited  in  the 
low  frequencies  by  the  distance  the  mass  car.  travel,  and  in  the  high 
frequencies  by  the  maximum  current  the  motor  coil  can  handle.  The 
maximum  peak  force  for  the  actuator  can  be  calculated  using  equation  3. 


-  m  w  d 
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wliere  m^^  is  the  mass  of  the  proof  mass,  w  is  the  frequency  of  the 
signal  driving  the  proof  mass,  and  d  is  the  maximum  proof  mass 
di splacement . 

Using  the  maximum  displacement  of  ±0.A2  inches  (0.16  inches  of  the 
nominal  stroke  is  lost  to  the  rubber  grommets  which  buffer  the  proof 
i7i<i.s.s  from  the  mounting  brackets).  This  results  in  a  maximum  actuator 
output  of  0.15  Ibf  for  the  fundamental  z-axis  bending  mode  (1.33  Hz), 
and  0.26  Ibf  for  the  fundamental  y-axis  bending  mode  (1.75  Hz).  The 
maximum  output  for  the  first  torsion  mode  is  13.3  Ibf,  far  beyond  the 
motor  coil  capability  of  4.2  Ibf.  To  verify  this  capability,  the 
actuator  was  driven  near  the  modal  frequency  with  the  base  fixed.  The 
input  amplitude  was  set  so  that  the  proof  mass  used  the  entire  range  of 
travel  without  hitting  the  end  brackets.  Table  2.1  shows  the  maximum 


force  outputs  and  efficiency  measured  for  each  actuator. 
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Actuator 

Frequency 

(Hz) 

1.25 

1.75 

1.25 

1.75 

1.25 

1.75 

1.25 

1.75 

Eff ici 

ency 

85. 

7% 

82. 

3% 

82. 

7% 

85. 

2% 

76. 

7% 

93. 

1% 

81. 

2% 

81. 

5% 

Table  2.1.  Actuator  Maximum  Force  Outputs 


Measured 
Maximum  (Ibf) 

0.114 

0.214 

0.110 

0.224 

0.102 

0.242 

0.108 

0.212 


Applying  the  efficiency  for  actuators  #2  and  #4  to  the  z-axis 
fundamental  bending  frequency  and  summing  the  maximum  forces  results  in 
a  maximum  control  force  of  0.246  Ibf.  Similarly,  actuators  #1  and  #3 
should  provide  0.462  Ibf  for  the  fundamental  y-axis  bending  frequency 
To  determine  the  ability  of  the  control  circuit  to  maintain  the 
inertial  position  of  the  proof  mass  when  the  actuator  is  uncommanded, 
the  the  slip  fable  was  driven  by  a  shaker  while  the  compensation  circuit 
was  operated  without  an  input  signal.  This  test  configuration  is  shown 
in  Figure  2.26.  The  acceleration  of  the  proof  mass  should  be  zero  if  it 
is  perfectly  floating.  If  not,  it  will  be  applying  a  force  to  the  beam, 
the  effect  of  which  is  determined  by  the  phase  of  the  response.  A  90° 
phase  lead  could  drive  the  beam  unstable  if  the  force  is  large  enough  to 
overcome  the  inherent  damping  in  the  beam.  An  out  of  phase  signal  would 
act  as  a  resonance,  leaving  the  structure  unaffected.  A  90°  phase  lag 
acts  as  a  viscous  damper,  and  an  in  phase  response  would  effectively  add 
inertia  to  the  system,  as  if  the  proof  mass  were  fixed  to  the  base. 
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The  base  motion  transfer  functions  measure  the  accleration  response 
of  the  proof  mass  to  the  force  applied  to  the  slip  table  by  the  shaker. A 
force  gauge  was  mounted  between  the  slip  table  and  shaker  armature  to 
determine  the  force  applied  to  the  base  and  a  sine  sweep  was  used  to 
drive  the  shaker  over  the  desired  frequency  range.  The  frequency 
response  functions  measured  for  each  actuator  are  shown  in  Figures  2.27 
through  2.30.  A  force  ratio  can  be  derived  from  the  magnitude  response 
by  adding  -6.1  dB.  Note  the  phase  of  the  response  is  approximately 
-90°,  which  indicates  the  uncommanded  actuators  should  act  essentially 
as  viscous  dampers  on  the  structure. 

To  estimate  the  damping  the  uncommanded  actuators  will  impart  to 
the  fundamental  bending  modes,  the  slip  table  was  driven  at  the  modal 
frequency  and  the  proof  mass  acceleration  measured.  Several  different 
force  levels  were  applied  to  the  slip  table  to  determine  the  variation 
of  proof  mass  acceleration  with  base  displacement  and  velocity.  The 
peak  base  velocities  were  calculated  using  equation  4,  and  peak 
displacements  were  derived  by  dividing  the  velocities  by  the  natural 
frequency . 


Base  Velocity 


386.4  in/sec 


where  is  the  force  driving  the  slip  table  base,  is  the  weight  of 

slip  table  base,  and  w  is  the  modal  frequency  in  rad/sec.  Results  of 
these  tests  for  the  fundamental  bending  frequencies  are  in  Tables  2.2 
through  2.5. 

The  force  out  of  the  uncommanded  actuators  is  not  a  constant  with 
respect  to  velocity,  which  would  have  represented  a  standard  viscous 
damper.  However,  a  first  order  approximation  can  be  made  to  the  damping 
by  assuming  an  average  value  for  each  actuator,  adding  the  averages  for 
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Figure  1 .11 .  Actuator  #1  Base  Motion  Transfer  Function 


Base  Weight 

F 

IN 

V 

B 

Dg 

(Ibf) 

(Ibf) 

(in/sec) 

(in) 

0.0841 

0.793 

0.101 

5.223 

0.1413 

1.331 

0.169 

0.2213 

2.085 

0.265 

0.2399 

1.081 

0.138 

10.921 

0.3589 

1.616 

0.206 

0.3936 

1.773 

0.226 

0.1035 

0.696 

0.063 

5.223 

0.2188 

1.472 

0.134 

0.3458 

2.387 

0.217 

0.3631 

1.167 

0.106 

10.921 

0.5188 

1.669 

0.152 

0.5957 

1.917 

0.174 

0.0204 


0.0302 


0.0427 


0.0255 


0.0347 


0.0368 


0.0174 


0.0277 


0.0403 


0.0237 


0.0299 


0.0332 


0.0257 


0.0227 


0.0205 


0.0238 


0.0215 


0.0207 


0.0250 


0.0187 


0.0169 


0.0203 


0.0179 


0.0173 


Peak  force  driving  the  base 
Peak  base  velocity 
Peak  base  displacement 

Force  applied  by  the  proof  mass  to  the  base 


Base  Weight 
(Ibf) 


5.223 


10.921 


5.223 


10.921 


F 

V 

D 

F 

IN 

B 

B 

PM 

(Ibf) 

(in/sec) 

(in) 

(Ibf) 

0.0462 

0.436 

0.055 

0.0153 

0.0708 

0.667 

0.085 

0.0195 

0.1585 

1.493 

0.190 

0.0324 

0.1622 

0.731 

0.093 

0.0221 

0.2754 

1.241 

0.158 

0.0292 

0.3350 

1.509 

0.192 

0.0347 

0.0955 

0.643 

0.058 

0.0197 

0.1429 

0.961 

0.087 

0.0240 

0.2958 

1.986 

0,181 

0.0347 

0.1641 

0.528 

0.048 

0.0148 

0.3589 

1.155 

0.105 

0.0257 

0.5129 

1.650 

0.150 

0.0292 

(jj 

n 

Base  Weight 

F 

IN 

V 

B 

Dg 

F 

PM 

HSH 

(Hz) 

(Ibf) 

(Ibf) 

(in/sec) 

(in) 

(Ibf) 

0.0412 

0.388 

0.049 

0.0114 

0.0295 

5.223 

0.0989 

0.931 

0.119 

0.0206 

0.0221 

1.25 

0.1445 

1.362 

0.173 

0.0281 

0.0206 

0.0617 

0.278 

0.035 

0.0134 

0.0484 

10.921 

0.1059 

0.477 

0.061 

0.0149 

0.0312 

0.3467 

1.562 

0.199 

0.0294 

0.0188 

0.0484 

0.326 

0.030 

0.0103 

0.0317 

5.223 

0.1365 

0.918 

0.083 

0.0175 

0.0191 

1.75 

0.2344 

1.577 

0.143 

0.0256 

0.0179 

0.0794 

0.256 

0.023 

0.0103 

0.0404 

10.921 

0.3467 

1.116 

0.101 

0.0197 

0.0176 

0.4955 

1.594 

0.145 

0.0242 

0.0152 

F  ■  Peak  force  driving  the  base 

V  ■  Peak  base  velocity 

D  ■  Peak  base  displacement 

F  ■  Force  applied  by  the  proof  mass  to  the  base 

PM 
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Table  2.5.  Base  Motion  Actuator  Response  -  Actuator  #4 


Base  Weight 

F 

I  N 

V 

B 

Dg 

F 

PM 

(Ibf) 

(Ibf) 

(in/sec) 

(in) 

(Ibf) 

0.0531 

0.500 

0.064 

0.0149 

5.223 

0.1161 

1.094 

0.139 

0.0253 

0.1862 

1.754 

0.223 

0.0365 

0.0933 

0.420 

0.054 

0.0141 

10.921 

0.1622 

0.731 

0.093 

0.0208 

0.3020 

1.360 

0.173 

0.0315 

0.0617 

0.415 

0.038 

0.0136 

5.223 

0.1718 

1.156 

0.105 

0.0247 

0.2541 

1.710 

0.155 

0.0322 

0.1334 

0.429 

0.039 

0.0134 

10.921 

0.2661 

0.856 

0.078 

0.0205 

0.3631 

1.168 

0.106 

0.0250 

0.0298 


0.0231 


0.0208 


0.0334 


0.0285 


0.0231 


0.0327 


0.0214 


0.0188 


0.0313 


0.0240 


0.0214 


Peak  force  driving  the  base 
Peak  base  velocity 
Peak  base  displacement 

Force  applied  by  the  proof  mass  to  the  base 


the  actuator  pair  controlling  the  axis  being  analyzed,  and  applying  the 
resulting  coefficient  as  a  damping  element  between  the  end  node  of  the 
structure  and  a  fixed  node.  The  element  is  oriented  in  the  plane  of  the 
axis  being  evaluated.  This  procedure  was  implemented  using  the  finite 
element  model  for  the  beam  with  floating  proof  masses  configuration  (see 
Chapter  3,  Table  3.4). 

The  finite  element  code  forms  the  system  of  equations  for  the  model 
in  the  form  of  (6:2-26,2-27) 

M  x(f)  +  C  x(f)  +  K  x(f)  -  F(f)  (5) 


where  x(f)  are  the  nodal  displacements,  F(f)  is  a  periodic  forcing 
function  applied  at  given  nodes,  and  f  is  a  discrete  frequency  value. 
This  system  is  assumed  to  be  in  a  steady  state  condition  at  the  given 
frequency,  thus  the  velocities  and  displacements  are 


X  (f)  -  -(2wf)^  x 

X  (f)  -  i(2jrf)  X 


(6) 


where  i  =  (-1)^^^. 


The  system  of  equations  now  becomes 


which  is  solved  for  each  designated  frequency  to  generate  a  frequency 
response  function.  The  damping  can  be  estimated  from  the  frequency 
response  data  using  the  half  power  point  method  (7:96)  which  has  the 


where  and  are  the  frequencies  at  which  the  response  is  3  dB  below 
the  resonant  peak,  w  . 

n 

Actuators  #2  and  #4  have  been  designated  to  control  the  z-axis  of 
the  beam,  and  #1  and  #3  will  control  the  y-axis  (see  Figure  5.3). 
Averaging  the  damping  coefficients  measured  from  the  base  motion  data 
and  summing  the  coefficients  for  the  two  actuators  results  in  a  z-axis 
damping  coefficient  of  0.0495  Ibf/in/sec,  and  a  y-axis  coefficient  of 
0.0431  Ibf/in/sec.  These  coefficients  were  applied  to  the  appropriate 
finite  element  beam  models  and  the  structural  frequency  response 
calculated  for  each  axis.  Results  are  shown  in  Figures  2.31  and  2.32 
for  the  z  and  y  axes  respectively.  The  damping  ratios  predicted  by  thi 
model  for  the  fundamental  bending  modes  are  in  Table  2.6. 

Table  2.6.  Predicted  Modal  Damping  Due  to  Actuators 


Axis  Predicted  Damping  Coefficient 


While  this  is  substantially  more  damping  than  would  normally  be 
encountered  in  an  uncontrolled  structure,  there  will  still  be  enough 
margin  between  the  residual  damping  and  the  maximvim  actuator  capability 
to  demonstrate  the  controller  is  affecting  the  structure. 

Finally,  since  the  actuators  will  operate  in  orthogonal  pairs,  the 
will  experience  side  force  loading  as  the  orthogonal  axis  vibrates.  To 
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ensure  this  loading  does  not  adversely  affect  the  actuator  operation, 
actuator  transfer  functions  were  measured  while  driving  the  slip  table 
with  the  actuator  mounted  perpendicular  to  the  direction  of  slip  table 
motion.  Tests  were  conducted  by  using  the  shaker  to  drive  the  slip 
table  at  two  different  amplitudes  for  each  of  the  first  five  modal 
frequencies.  A  sine  sweep  was  used  to  generate  an  actuator  transfer 
function  over  the  control  bandwidth  while  the  slip  table  was  in  motion. 
Variation  of  the  actuator  response  with  off  axis  force  and  frequency  was 
minimal  and  no  trends  were  observed  in  the  response  variation.  The 
standard  deviation  from  the  nominal  transfer  function  of  all  the  tests 
appears  in  Figure  2.33,  and  the  variation  seems  to  be  mainly  measurement 
error.  The  stability  of  the  response  under  these  conditions  indicates 
that  the  actuators  should  not  be  affected  by  off-axis  forces. 


Comparison  of  Response  with  Design  Goals 

Table  2.7  lists  the  final  transfer  function  for  each  actuator  at 

each  modal  frequency.  The  response  still  has  some  magnitude  and  phase 

variation,  however  it  was  decided  the  magnitude  variation  could  be 

accounted  for  in  the  state  space  model  and  the  phase  variation  would  be 

ignored.  While  ignoring  20°  of  phase  at  the  fundamental  bending 

frequency  may  seem  excessive,  this  was  the  maximum  observed  for  all 

input  amplitude  levels.  A  problem  which  remains  unresolved  is  the  phase 

response  variation  with  command  signal  amplitude.  The  phase  was 

o 

observed  to  increase  as  much  as  15  when  the  input  amplitude  was  reduced 


by  half.  The  data  presented  in  the  '  ible  2.7  is  for  an  input  amplitude 
that  gives  the  maximum  force  output  from  the  actuator  near  the  z-axis 
fundamental  bending  frequency.  Any  reduction  in  force  command  will 
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Figure  2.33.  Standard  Deviation  of  the  Variation  in  Actuator  Frequency  Response  with  Side  Force  Loading 


Table  2.7.  Final  Actuator  Transfer  Functions 


Actuator 

Frequency  (Hz) 

1.25 

1.75 

13.0 

19.0 

25.25 

*#1 

* 

Mag 

0.5828 

0.4904 

0.4126 

0.4032 

0.4032 

Phase* 

-19.3 

-16.1 

4.9 

6.6 

8.5 

#2 

Mag 

0.5893 

0.4959 

0.3939 

0.3939 

0.3984 

Phase 

-20.9 

-16.2 

4.2 

7.3 

10.5 

#3 

Mag 

0.5800 

0.4880 

0.4153 

0.4106 

0.4059 

Phase 

-18.7 

-15.2 

5.8 

7.1 

8.2 

#4 

Mag 

0.5858 

0.4929 

0.3961 

0.3961 

0.4007 

Phase 

-21.7 

-17.2 

4.6 

8.3 

10.9 

*  Magnitude  in  Ibf/V  #  Phase  in  degrees 


increase  the  phase  response,  thus  moving  it  closer  to  zero.  A 
discussion  and  characterization  of  this  problem  is  in  Appendix  B. 

The  actuators  do  add  some  damping  to  the  structure,  thus  the  proof 
masses  are  not  perfectly  floating,  however  the  level  is  low  enough  that 
showing  significance  between  residual  and  commanded  damping  levels 
should  not  present  a  problem.  The  greater  concern  is  the  point  at  which 
the  proof  masses  come  to  rest  and  suddenly  contribute  to  the  system 
inertia.  Uhile  this  amplitude  level  is  small,  and  may  be  reduced  even 
more  when  the  actuators  are  commanded,  it  makes  testing  with  the  modal 
analyzer  difficult,  since  the  response  appears  as  two  separate  modes 
which  are  present  and  different  times. 

The  proof  masses  do  maintain  a  centered  position  when  commanded, 
though  some  change  in  center  still  was  observed  with  input  amplitude. 
With  the  sensitivity  to  bias  shifts  that  these  actuators  exhibit,  the 
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shift  could  easily  be  in  the  frequency  generator  of  the  Ono  Sokki .  A 
3  mV  bias  shift  in  a  2  V  signal  is  a  very  small  percentage  error  and 
would  be  difficult  to  detect. 
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No  specific  measures  were  taken  to  limit  the  proof  mass  travel. 

The  structural  deflections  will  be  limited  to  prevent  the  actuators  from 
saturating. 


Unresolved  Problems,  Conclusions,  and  Recommendations 


The  variation  in  the  final  actuator  transfer  function  will  induce 
some  errors  if  it  is  ignored  in  the  overall  system  model.  To  minimize 
this  error,  the  response  magnitude  can  be  accounted  for  in  the  state 
space  model,  however  the  error  induced  by  the  phase  variation  is 
probably  best  ignored  since  it  will  vary  with  input  amplitude  and  would 
significantly  increase  the  complexity  of  the  model.  Future  work  to 
eliminate  this  problem  is  recommended  to  increase  the  accuaracy  of  the 
actuator  response  and  the  overall  experiment.  Another  problem  which 
remains  unresolved  is  noise  in  the  system,  especially  at  60  Hz. 

Filtering  of  noise  this  close  to  the  control  bandwidth  is  difficult 
since  the  filter  dynamics  will  affect  the  response  of  the  system. 
Filtering  the  power  inputs  to  the  compensation  circuit  may  provide  some 
isolation  and  reduce  the  noise  amplitude. 

Several  concerns  about  linear  proof  mass  actuators  as  viable 
control  mechanisms  for  this  experiment  and  for  large  space  structures  in 
general  have  arisen  from  the  development  for  this  experiment.  First  is 
the  limited  force  output  at  very  low  frequencies.  Since  the  force 
output  is  derived  from  accelerating  a  mass,  large  low  frequencv  force 
outputs  require  either  a  very  long  travel  for  the  mass  or  a  large  mass. 


i 
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Both  approaches  have  problems  when  considered  for  space  structure 
control.  If  the  proof  mass  is  given  a  long  travel,  linearity  of  the 
motor  output  becomes  a  problem.  Perhaps  this  nonlinearity  could  be 
modeled  and  accounted  for,  but  the  cost  would  be  a  significant  increase 
in  the  computational  requirements  on  the  control  system.  Large  masses 
present  two  disadvantages.  First,  unless  the  moving  mass  is  designed  to 
consist  at  least  partially  of  hardware  that  would  normally  be  a  part  of 
the  vehicle  the  increased  mass  is  merely  dead  weight  which  must  be 
boosted  into  orbit.  Second,  and  •"ore  significant  to  the  control 
probleiT.,  is  adding  large  masses  to  a  light,  flexible  structure  tends  to 
create  nodal  points  at  or  very  close  to  the  actuators.  The 
significantly  reduced  modal  amplitude  at  the  control  mechanism  requires 
the  actuator  to  deliver  more  force  to  maintain  the  same  level  of 
control.  This  problem  is  not  as  significant  if  the  proof  mass  is 
controlled  such  that  it  is  not  contributing  to  the  structural  mass  (as 
is  the  case  for  this  experiment),  but  if  the  actuators  are  mounted  on 
the  structure  in  colocated,  orthogonal  pairs,  the  orthogonal  actuator 
does  create  a  node.  Future  experiments  and  actual  structures  which  use 
linear  proof  mass  actuators  should  avoid  colocating  the  actuators  if 
their  mass  is  large  in  relation  to  the  structure  being  controlled. 
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III.  Structure 
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Description 

The  structure  to  be  controlled  in  this  experiment  is  a  cantilevered 
beam  with  rectangular  cross  section.  A  circular  plate  is  mounted  on  the 
free  end  of  the  beam  to  provide  a  surface  to  mount  the  actuators  (see 
Figure  3.1).  The  configuration  and  dimensions  are  shown  in  Figure  3.2 
and  Table  3.1  lists  the  physical  charcteristic  of  the  structure.  This 
configuration  was  selected  by  the  Flight  Dynamics  Laboratory  to  exhibit 
the  large  space  structure  charcteristics  of  low  frequency,  closely 
spaced  vibration  modes.  It  also  allowed  for  ungrounded  sensing  and 
actuation  which,  with  vertical  suspension  of  the  beam,  provides  an 
approximation  to  the  zero  gravity  free  vibration  environment  of  a  large 
space  structure. 

.*iodelling 

There  are  many  modelling  methods  available  to  find  the  natural 
vibration  frequencies  and  mode  shapes  for  structures.  The  cantilevered 
beam  with  a  tip  mass  is  a  simple  enough  structure  to  permit  analytical 
solution  of  the  free  vibration  equations  of  motion.  Other  techniques 
which  form  approximations  to  this  exact  solution  include,  but  are  not 
limited  to,  Rayliegh-Ritz ,  Assumed  Modes,  Collocation,  Galerkin,  and 
Finite  Element  methods.  Because  the  Finite  Element  Method  is  highly 
automated  and  easily  applied  it  was  used  as  the  primary  modelling  tool 
for  this  structure.  To  verify  the  validity  of  the  solution,  the 
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Figure  3.2.  Structure  Configuration 


Property  Description 


Value 


Units 


Beam  Length  (L) 

70.75 

in 

Y  Cross-Section  Width  (a) 

1.01 

in 

Z  Cross-Section  Width  (b) 

0.758 

in 

Cross-Section  Area  (A) 

0.7656 

.  2 
in 

Young's  Modulus  (E) 

10.8x10® 

psi 

Shear  Modulus  (G) 

4.1x10® 

psi 

Beam  Density  (p  ) 

B 

2.591x10''' 

Ibf  •  sec  ®/in'' 

Beam  Mass  (M  ) 

B 

1.403x10'^ 

Ibf'sec  /in 

Y  Moment  of  Inertia  (I  ) 

y 

3.667x10'^ 

.  A 

in 

Z  Moment  of  Inertia  (I  ) 

Z 

6,508x10'^ 

.  A 

in 

Torsional  Moment  of  Inertia  (K) 

7.913x10'^ 

.  A 

in 

Polar  Mass  Moment  of  Inetia  (I  ) 

m  X 

1.865x10'® 

2 

Ibf'sec  /in 

Plate  Diameter  (d) 

12.0 

in 

Plate  Thickness  (t) 

1.0 

in 

Plate  Mass  (M^)  (measured) 

2.847x10'^ 

2 

Ibf'sec  /in 

X  Mass  Moment  of  Inertia  (I  ) 

D  m  X 

0.5125 

2 

Ibf'sec  'in 

Y-Z  Mass  Moment  of  Inertia  (I  ) 

Dm 


Ibf • sec  • in 


0.2562 


equations  of  motion  were  solved  for  a  cantilevered  beam  with  tip  mass 
and  compared  with  the  finite  element  results.  The  finite  element  model 
was  implemented  using  the  MSC*PAL  software  package  availible  at  the 
Flight  Dynamics  Laboratory.  MSC*PAL  is  a  finite  element  code  which  runs 
on  an  IBM-PC  or  compatible  micro-computer. 

The  structure  was  modelled  in  two  steps.  First,  the  finite  element 
and  exact  solutions  were  formulated  and  compared  for  the  beam  without 
the  actuators  mounted.  The  finite  element  solution  was  modified  by 
adding  elements  to  the  model  until  the  solution  converged  and  matched 
the  solution  to  the  equations  of  motion.  Using  this  finite  element 
model  the  mass  model  for  each  actuator  was  added  to  determine  the 
natural  frequencies  and  mode  shapes  for  the  controlled  configuration. 
Since  an  ideal  actuator  will  have  a  floating  proof  mass,  this  mass  will 
not  contribute  to  the  tip  mass  in  the  actuator  control  axis.  Thus,  the 
final  model  is  actually  a  composite  of  three  different  models,  one  for 
the  tip  mass  configuration  of  each  bending  direction  and  torsion. 

The  classical  solution  for  bending  vibration  of  beams  can  be  found 
in  many  textbooks.  The  equations  of  motion  are  formulated  using 
Hamilton's  Principal,  and  after  applying  separation  of  variables,  the 
resulting  boundary  value  problem  can  be  solved  for  the  natural 
frequencies  and  mode  shape  of  the  structure  (9:161)  In  general,  shear 
and  rotatory  inertia  effects  are  ignored  for  beams  with  cross-sectional 
areas  that  are  small  compared  with  the  length,  thus  the  equations  of 
motion  for  bending  vibration  of  this  structure  are 
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Young's  Modulus  for  the  beam 

Moment  of  inertia  of  the  beam  in  the  bending  axis 
Displacement  function  of  the  beam 
Distance  along  the  beam 
Mass  per  unit  length  of  the  beam 


Separating  variables  and  applying  boundary  conditions  results  in  a 
characteristic  equation  that  is  solved  numerically  for  it  roots,  which 
are  the  natural  vibrations  frequencies  of  the  beam.  Through  a  similar 
procedure  the  torsional  equation  of  motion  is  formulated  (9:156)  and  is 
found  to  be 
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where  GJ  ■  Torsional  rigidity  of  the  beam 

I  ■  Polar  mass  moment  of  inertia  of  the  beam 
6  s  Angular  displacement  of  the  beam 

Again,  separating  variables  and  applying  boundary  conditions  results  in 
a  charcteristic  equation  which  can  be  solved  numerically  for  the  natural 
frequencies . 

The  Finite  Element  Method  formulates  the  problem  by  breaking  the 
structure  into  discrete  elements  and  assuming  a  form  of  the  displacement 
function  over  each  element.  The  elements  are  connected  at  nodal 
points  and  thus  form  a  piecewise  approximation  to  the  structure  (10:2). 
Each  element  has  a  mass  and  stiffness  matix  associated  with  it,  which  is 
defined  by  the  beam  properties  and  the  specific  displacement  function 
assumed  for  the  element.  Combining  these  elemental  matrices  results  in 
a  system  of  equations  of  the  form 
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Mass  matrix 

Stiffness  matrix 

Nodal  displacement  vector 

Vector  of  forces  applied  at  the  nodes 


(13) 


The  homogeneous  solution  to  this  coupled  system  of  equations  represents 
the  free  vibration  of  the  structure  from  which  the  natural  frequencies 
and  mode  shapes  are  found  by  solving  for  the  eigenvalues  and 


eigenvectors.  The  accuracy  of  this  method  when  applied  to  a  simple  beam 
depends  primarily  on  the  number  of  elements  chosen  to  model  the 
structure.  For  this  case,  a  10  element  model  was  found  to  match  the 
exact  solution  for  the  lower  modes  to  reasonable  accuracy.  The  higher 
modes  will  show  increasingly  more  error  as  the  mode  number  approaches 
the  number  of  elements  in  the  model.  Table  3.2  compares  the  natural 
frequencies  for  the  structure  that  were  derived  from  the  two  solutions. 

Using  the  finite  element  model  for  the  basic  beam,  the  actuator 
mass  models  were  added  to  predict  the  natural  frequencies  and  mode 
shapes  for  the  structure  in  the  control  configuration.  The  actuator 
mass  models  were  developed  by  breaking  them  into  six  components  which 
could  be  approximated  by  simple  figures.  The  mass,  center  of  gravity, 
and  mass  moments  of  inertia  for  each  of  these  figures  were  then  used  in 
the  finite  element  model.  The  MSC-PAL  software  supports  adding  a  lumped 
mass  at  some  point  on  the  structure  and  assigning  the  effects  of  that 
mass  to  a  specific  node  point.  All  actuator  mass  elements  were  assigned 
to  the  end  node  of  the  beam  (#11),  along  with  the  plate.  The  software 
automatically  uses  the  parallel  axis  therom  to  calculate  the  moments  of 
inertia  for  each  mass  element  about  the  assigned  node.  The  detailed 
model  for  each  actuator  can  be  found  in  Appendix  C.  To  approximate  a 
floating  proof  mass,  the  mass  elements  representing  the  proof  masses 
were  omitted  for  the  actuators  controlling  the  modelled  axis.  For 
torsion,  all  proof  masses  were  omitted  since  for  small  displacements 
they  should  not  contribute  to  the  structure  mass  properties.  Natural 
frequencies  from  the  beam/actuator  finite  element  model  are  shown  in 


Table  3.3. 


Mode  Axis 


Modal  Frequencies  (Hz) 

FEM 

Exact 

Finite  Element 

error  (%) 

1.611 

1.608 

0.19 

2.146 

2.136 

0.47 

14.401 

14.401 

0.00 

20.968 

20.749 

1.04 

27.938 

27.642 

1.06 

56.974 

57.115 

0.25 

75.911 

76.092 

0.24 

103.382 

104.656 

1.2  3 

137.744 

139.414 

1.21 

180.088 

181.375 

0.71 

Modal  Frequency 


1.325 


1.773 


12.796 


19.431 

25.883 


47.896 


63.817 

86.451 


Modal  Tests 


The  most  accurate  method  for  determining  the  frequencies  and  mode 
shapes  of  a  structure  is  almost  always  modal  testing.  While  models  give 
insight  into  the  phenomena  being  observed  in  a  system,  they  are  still 
only  mathematical  representations  of  reality.  Unfortunately,  modal 
testing  may  be  impractical  for  many  large  space  structures.  Due  to 
their  size  and  flexibility,  they  may  not  even  be  able  to  be  assembled 
until  they  are  on  orbit.  Even  if  they  can  be  assembled,  simulating  a 
zero  gravity  environment  for  a  large  structure  may  prove  extremely 
difficult.  WTiile  actual  space  structures  do  present  these  problems, 
maintaining  the  accuracy  of  the  experiment  is  a  paramount  consideration 
and  the  structura  model  for  the  control  design  will  use  modal  test  data 
whenever  possible. 

Modal  tests  were  conducted  using  the  impact  method  which  excites 
the  structure  by  tapping  it  with  a  force  hammer.  The  hammer  has  a 
piezoelectric  crystal  in  the  tip  which  registers  a  voltage  signal  when 
stressed.  The  structural  response  to  this  force  is  measured  using 
accelerometers  located  at  various  points  on  the  structure.  The  force 
and  accelerometer  signals  are  recorded  by  a  GenRad  Modal  Analyzer,  which 
Fast  Fourier  Transforms  the  time  response  data  and  generates  a  frequency 
response  function  for  the  structure.  Natural  modes  of  vibration  for  the 
structure  show  up  as  peaks  in  this  frequency  response.  Mode  shapes  can 
be  determined  by  measuring  the  response  at  several  different  points 
along  the  beam.  The  amplitude  and  direction  of  the  mode  at  each  point 
is  defined  by  the  amplitude  and  phase  of  the  frequency  response  mea.sur*  d 
the  re . 

The  first  modal  test  was  conducted  on  the  basic  beam  configuration 
without  the  actuators  Results  of  this  tes’  and  the  difference  between 


the  measured  response  and  finite  element  model  are  shown  in  Table  3.4. 
While  the  fundamental  frequencies  match  fairly  well  with  the  model,  the 
errors  grow  unexpectedly  large  for  the  second  bending  modes  and  beyond. 
The  modelling  of  a  cantilevered  beam  is  a  straight  forward  process  and 
should  generate  accurate  predictions  well  beyond  the  second  mode.  To 
evaluate  this  problem,  Mr.  Bob  Gordon  of  AFWAL/FIBG  developed  a  series 
of  models  which  varied  both  the  beam  parameters  and  the  boundary 


Table  3.4.  Measured  vs  Finite  Element  Modal  Frequencies  -  No  Actuators 


y. 


Mode 

Axis 

Modal  Frequencies  (Hz) 

FEM 

error  (%) 

Measured 

Finite  Element 

1 

Z 

1.617 

1.608 

0.56 

2 

Y 

2  . 129 

2.136 

0.33 

3 

T 

14.436 

14.401 

0.24 

4 

Z 

20.160 

20.749 

2.92 

5 

Y 

26.847 

27.642 

2.96 

Bi 

Z 

54.828 

57.115 

4.17 

■ 

Y 

72.985 

76.092 

4 , 26 

■ 

Z 

98.932 

104.656 

5.79 

conditions  of  the  model.  Results  of  this  effort  were  inconclusive. 
Whiilc’  a  more  accurate  prediction  of  the  high  frequency  modes  could  be 
achieved,  it  was  always  at  the  cost  of  adding  error  to  the  fundamental 
modes.  Neither  changing  beam  parameters,  nor  changing  the  boundary 
condition  from  a  perfect  cantilever  to  a  pinned  joint  with  a  rotational 
spring  had  the  'effect  of  maintaining  the  fundamental  frequencies  while 
lowering  the  higher  mode  predicted  frequencies.  Mr.  Gordon  also  found 
that  a  model  of  the  entire  support  frame  and  beam  resulted  in  a  verv 
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ill-conditioned  system  of  equations  which  induced  substantial  numerical 
errors  into  the  solution.  The  ill-conditioning  results  from  the  support 
plate  to  which  the  beam  is  attached  (see  Figure  3.1).  This  plate  is  the 
most  massive  and  stiffest  element  in  the  frame/beam  structure,  while  the 
least  massive  and  stiff  element  is  the  attached  beam.  While  the 
difference  between  the  model  and  modal  test  data  probably  lies  in  some 
coupling  between  beam  and  support  structure,  it  could  not  be 
demonstrated  in  the  model. 

The  model  was  considered  acceptable  for  the  first  five  modes, 
consisting  of  the  first  two  bending  modes  in  each  direction  and  the 
first  torsion  mode,  since  the  prediction  of  the  fundamental  modes  is 
fairly  accurate  and  the  errors  in  the  second  bending  modes  were  less 
than  3%.  While  the  remaining  modes  exhibited  larger  errors,  they  are 
outside  the  control  bandwidth  for  the  experiment  and  as  such  will  not  be 
incorporated  in  the  controller. 

A  second  modal  test  was  performed  with  the  actuators  mounted  on  the 
beam.  This  test  was  performed  prior  to  completion  of  the  actuator 
compensators;  therefore,  to  simulate  the  ideal  actuator,  the  proof 
masses  were  removed  from  the  actuators  controlling  the  axis  under  test. 
For  torsion,  all  proof  masses  were  removed.  Results  of  this  modal  test 
are  in  Table  3.5.  The  data  shows  the  finite  element  model  for  this 
configuration  to  have  errors  similar  to  the  basic  beam  tests  for  the 
first  five  modes.  Errors  are  more  substantial  in  the  higher  modes. 

Several  problems  were  encountered  in  conducting  the  modal  tests 
which  relate  to  low  frequency  testing  of  lightly  damped  structures. 
First,  the  lighter  the  damping,  the  more  difficult  it  is  to  accurately 
estimate  a  damping  coefficient  for  the  mode.  The  software  in  the  GenRad 
system  estimates  damping  by  fitting  a  complex  exponential  curve  to  the 
frequency  response  data.  To  accurately  fit  the  data  for  a  lightly 


damped  structure  requires  very  small  frequency  steps  in  the  transformed 
data.  If  the  frequency  step  is  too  large,  the  true  resonant  peak  will 
be  aliased  out  of  the  data.  In  addition,  the  peak  that  is  captured  may 
be  so  sharp  that  a  single  frequency  step  may  drop  below  the  half  power 


Table  3.5.  Measured  vs  Finite  Element  Modal  Frequencies  -  With  Actuators 


Mode 

Axis 

Modal  Frequencies  (Hz) 

FEM 

error  (%) 

Measured 

Finite  Element 

1 

Z 

1.3A6 

1.325 

1.56 

2 

Y 

1.761 

1.773 

0.68 

3 

T 

12.963 

12.796 

1.29 

A 

Z 

18.877 

19.431 

2.93 

5 

Y 

25.141 

25.883 

2.95 

6 

Z 

44.567 

47.896 

7.47 

7 

Y 

59.861 

63.817 

6.81 

8 

Z 

86.451 

94.875 

9.74 

points.  Without  these  data  points,  the  accuracy  of  the  damping  estimate 
becomes  very  suspect.  This  situation  existed  for  both  of  the  modal 
tests.  Even  with  the  GenRad's  smallest  frequency  step  the  true  resonant 
peaks  and  half  power  points  were  rarely  captured.  Damping  ratio 
estimates  for  most  modes  were  below  0.0005  and  were  in  actuality 
probably  less. 

The  second  problem  this  caused  was  to  make  accurate  measurement  of 
the  mode  shapes  almost  impossible.  Missing  the  true  resonant  peak  of 
the  mode  corrupts  the  estimate  of  the  modal  amplitude.  In  addition,  the 
change  in  the  phase  of  the  response  occurs  so  quickly  that  one  frequency 
step  shifts  the  phase  180°,  effectively  missing  the  ±  90°  phasing  of  the 
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The  final  modal  test  was  accomplished  with  the  structure  and 
actuators  in  the  final  configuration  for  control  testing.  All  wiring 
for  sensors  and  actuators  was  attached  to  the  beam  and  the  compensated 
actuators  were  mounted  and  operating  without  command  inputs.  Data 
collection  and  analysis  was  again  accomplished  with  the  GenRad  system, 
with  the  exception  of  the  fundamental  bending  modes.  One  of  the 
problems  with  the  actuators  was  the  tendency  of  the  proof  masses  to 
stick  when  the  amplitude  of  the  oscillation  in  the  fundamental  modes 
became  small.  The  result  is  the  structure  effectively  has  two  closely 
spaced  modes,  one  when  the  proof  masses  are  floating  and  a  second  when 
they  become  fixed.  Because  of  this,  the  GenRad  proved  incapable  of 
providing  good  estimates  of  the  fundamental  bending  modes.  To  test 
these  modes,  a  sine  dwell  was  conducted  by  fixing  small  magnets  to  the 
base  plate  in  the  axis  to  be  tested  and  forcing  the  beam  at  a  single 
frequency  using  magnetic  induction  coils.  A  frequency  response  function 
was  generated  by  incrementally  increasing  the  frequency  of  the  forcing 
function,  allowing  the  motion  to  become  steady  state,  and  measuring  the 
magnitude  and  phase  of  the  beam  acceleration  with  respect  to  the  forcing 
signal.  The  resonant  peak  of  this  response  function  was  used  to  define 
the  modal  frequency  and  the  modal  damping  was  estimated  using  the  half 
power  points  of  resonant  peak  (see  equation  8).  The  mode  shape  was 
estimated  by  forcing  the  scruc*:ure  at  the  resonant  freqeuency  and 
measuring  the  acceleration  at  several  points  along  the  beam.  The 
stations  used  for  this  test  coincide  with  the  node  points  of  the  finite 
element  model,  which  allows  for  consistent  comparison  of  results  with 
the  finite  element  model.  This  data  and  the  data  for  the  higher  modes 
obtained  using  the  GenRad  can  be  found  in  Appendix  D. 

A  summary  of  the  natural  frequencies  is  presented  in  Table  36  and 
tbe  mode  shape.s,  normalized  with  respect  to  the  mass  matrix,  are  shown 
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Table  3.6.  Modal  Frequencies  and  Damping  for  Control  Configuration 


Mode 

Axis 

Modal  Frequencies  (Hz) 

FEM 

error  (%) 

Damping 

Measured 

Finite  Element 

1 

Z 

1.33 

1.329 

0.08 

0.0377 

2 

Y 

1.75 

1.775 

1.43 

0.0350 

3 

T 

12.580 

12.655 

0.60 

0.0127 

4 

Z 

18.094 

19.394 

7.21 

0.0029 

Y 

23.897 

25.839 

8.11 

0.0032 

■ 

Z 

43.454 

47.734 

9.86 

0.0022 

7 

Y 

57.624 

63.608 

10.39 

0.0014 

in  Table  3.7.  Since  the  generalized  mass  was  not  derived  from  the  modal 
test  data  the  value  calculated  in  the  finite  element  model  was  used  to 
appropriately  scale  the  measured  mode  shapes.  The  error  induced  should 
be  small  since  the  model  closely  matched  the  first  five  modes.  The 
actuators  have  added  enough  damping  to  the  structure  to  alleviate  the 
problem  previously  encountered  with  estimating  damping  and  mode  shapes. 
The-  finite  element  modal  was  found  to  give  a  slightly  better 
approximation  to  the  modal  test  data  if  the  inertia  of  the  proof  masses 
was  included  without  their  mass.  The  correction  mainly  affected  the 
fundamental  modes  and  slightly  improved  the  higher  modes. 

The  normalized  bending  mode  shapes  measured  for  the  first  two 
Ijerding  modes  in  each  axis  are  compared  with  the  finite  element  model 
mode  shapes  in  Figures  3,4  through  3.7.  The  mode  shape  for  the 
torsional  vibration  mode  was  not  measured  in  the  modal  test  because  the 


beam  was  not  instrumented  to  measure  torsion  at  anv  location  other  tnan 
the  base  pi  ate  Since  the  finite  element  fundamental  bending  modes 

showed  good  agreement  with  the  modal  tests,  and  the  torsion  frequency 
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Figure  3.4. 


2-Axis  Fundamental  Bending  Mode 
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Figure  3.6.  Z-Axis  Second  Bending  Mode 


Station  Number 


Figure  3.7.  Y-Axis  Second  Bending  Mode 


Figure  3.8.  First  Torsion  Mode  -  Finite  Element  Model  Only 


IV.  Sensors 


In  keeping  with  the  goal  of  constructing  an  experiment  with  as  much 
in  common  with  large  space  structures  as  possible,  the  sensors  selected 
for  use  in  controlling  the  beam  were  piezoelectric  accelerometers. 
Accelerometers  provide  an  inertial  sensing  capability  which  is  required 
in  space.  The  primary  disadvantages  of  using  accelerometers  are  that 
the  output  must  be  integrated  to  generate  a  velocity  measurement  for  the 
control  algorithm,  and  that  most  accelerometers  of  this  type  have  not 
been  designed  for  an  accurate  response  below  1  Hz.  While  integration 
does  tend  to  smooth  some  high  frequency  noise,  it  also  aggrevates  very 
low  frequency  drift  and  bias  signals.  An  example  of  the  drift  seen  in 
one  accelerometer  is  shown  in  Figure  4.1. 

The  specific  accelerometers  selected  for  sensors  on  the  beam  were 
manufactured  by  Kistler  Instrument  Corporation.  Specifications  for 
these  devices  are  listed  in  Appendix  E.  These  instruments  had  the  best 
low  frequency  response  characteristics  and  the  least  drift  of  the 
accelerometers  available  at  the  Flight  Dynamics  Laboratory  for  this 
experiment.  Even  so,  the  output  was  found  to  have  so  much  low  frequency 
drift  the  actuators  would  not  remain  centered  when  a  control  loop  was 
closed.  Because  the  complex  pole  pair  in  the  actuator  dynamics  has  been 
moved  below  1  Hz,  they  are  verv  responsive  to  these  low  freqvienrv 
sigr.al.s  As  was  noted  earlier,  a  b  mV  drift  in  the  input  signal  hat; 


hieen  observed  to  drive  the  proof  mass  to  the  end  of  it.s  stroke 


To  remove  this  drift  from  the  velocity  measurements  the 
accelerometer  output  was  passed  through  a  signal  conditioning  circuit 
prior  to  integration.  Ideal  signal  conditioning  for  this  experiment 
would  reject  all  signals  below  1  Hz  and  not  affect  the  rest  of  the 
frequency  spectrum.  Unfortunately,  filters  which  reduce  the  magnitude 
of  the  signal  affect  the  phase  for  a  decade  beyond  the  break  frequency. 
Therefore,  to  ensure  that  the  phase  of  the  signal  is  unaffected  at  the 
fundamental  bending  frequency  of  the  beam  requires  placing  the  break 
frequency  at  or  below  0.15  Hz.  The  problem  again  arises  that  noise  in 
the  system  is  so  close  in  frequency  to  the  control  bandwidth  that 
filtering  dynamics  contaminate  the  measurement  response  in  this 
bandwidth . 

The  approach  taken  to  filter  some  of  the  low  frequency  input  was  to 
split  the  signal  and  on  one  side,  low-pass  filter  it  to  remove 
frequencies  above  1  Hz,  invert  the  filtered  signal  and  then  sum  it  with 
the  original,  unfiltered  signal.  This  should  have  the  effect  of 
subtracting  out  the  drift.  The  filter  and  integrator  were  implemented 
using  the  circuit  in  Figure  4.2.  The  resulting  transfer  function  from 
acceleration  to  velocity  is 


>:(s) 

>:(s) 


3.5081  (S^  ^  +  7.814  ±  10.76j)(S^  ^  +  20.46  ±  6.628j) 
(S  +  11.31)^  (S  -f  0.6378) 


(14) 


The  predicted  frequency  response  for  this  circuit  is  in  Figure  4.3. 
This  filtering  has  reduced  the  response  magnitude  of  the  0.1  to  1  Hz 
frequency  range  when  compared  with  a  simple  integrator.  The  transfer 


function  for  an  integrator  has  a  15  dB  drop  in  magnitude  from  0.2  Hz  to 
1.3  Hz.  The  peak  magnitude  for  the  filter  occurs  at  0.4  Hz  and  is  only 
7  dB  above  the  1.3  Hz  response  magnitude.  While  it  would  be  desireable 
to  move  the  low  frequency  cutoff  even  higher,  the  phase  response  shows 
that  the  phase  has  just  settled  at  -90°,  so  moving  the  filter  higher 
would  begin  to  shift  the  phase  of  the  output  away  from  a  velocity 
signal.  The  frequency  response  of  the  actual  circuit  comes  very  close 
to  the  prediction,  as  can  be  seen  in  Figures  4.4  through  4.7.  As  a 
final  step  in  the  signal  conditioning  process  the  output  from  the 
integrator  was  passed  through  an  A/C  coupled  Intek  Variable  Gain 
Amplifier  to  remove  any  bias  that  may  still  be  present  in  the  signal. 

The  break  frequency  of  this  A/C  coupling  is  at  0.16  Hz,  which  accounts 
for  the  extra  phase  shift  seen  in  the  measured  frequency  response  verses 
prediction. 

Sensor  Calibration  and  Scale  Factors 

The  accelerometers  were  calibrated  using  a  1  g  peak  shaker  which 
operates  at  80  Hz.  Ideally  a  frequency  response  function  would  have 
been  generated  for  each  sensor  over  the  control  bandwidth. 

Unfortunately,  the  only  shakers  available  of  high  enough  quality 
response  to  be  used  as  a  calibration  instrument  had  a  low  frequency 
cutoff  of  10  Hz.  As  a  result,  each  sensor  was  calibrated  at  a  single 
frequency  and  the  manufacturer  specification  of  a  flat  response  function 
over  the  operating  range  of  the  sensor  was  relied  upon.  The  calibration 
measured  for  each  accelerometer  is  listed  in  Table  4.1. 
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Table  4.1.  Sensor  Calibration  Data 


Serial  # 

Sensor  » 

Sensitivity  (V/g) 

C82629 

1 

0.881 

C82753 

2 

0.881 

C82736 

3 

0.891 

C82702 

1 _ 

4 

0.871 

To  achieve  the  best  accuracy  possible  for  the  sensor  measurements 
the  response  functions  measured  for  the  signal  conditioning  and 
integrator  (Figures  4.4  to  4.7)  were  adjusted  to  match  the  predicted 
response  magnitude  at  1.3  Hz  by  adding  or  subtracting  the  required 
gain.  The  nominal  and  adjusted  gains  for  each  sensor  channel  are  listed 
in  Table  4.2.  Combining  these  gains  with  the  accelerometer 
sensitivities  results  in  the  sensor  scale  factors  shou-n  in  Table  4.3. 


Table  4.2.  Integration  and  Signal  Conditioning  Gain 


fel>] 


mV/ft/sec 

92.19 


mV/in/sec 

7.68 


2 

93.26 

1 .11 

3 

90.10 

7.51 

90.10 


7,51 
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System  Configuration  and  Gain  Verification 
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The  hardware  used  in  the  experiment  consists  of  the  beam  and  four 
proof  mass  actuators,  four  accelerometers  as  sensors,  several  stages  of 
signal  conditioning  on  inputs  and  outputs,  and  the  PC-1000  Systolic 
Array  Processor  and  its  host  computer,  a  Compaq  Portable  Computer.  Data 
collection  was  accomplished  using  the  Ono  Sokki  Freqeuncy  Spectrum 
Analyzer  with  an  HP-7874  Plotter,  and  the  GenRad  Modal  Analyzer  (see 
Figures  5.1  and  5.2). 

The  actuators  were  mounted  on  the  base  plate  in  orthogonal  pairs 
parallel  to  the  y  and  z  beam  axes  (see  Figure  5.3).  This  configuration 
provides  symmetric  forces  from  each  pair  for  bending  control  and 
asymmetric  forces  from  all  four  for  torsion  control.  The  actuator 
located  along  the  minus  z  direction  was  arbitrarily  selected  as  «1 ,  and 
remaining  actuators  were  numbered  from  there  in  a  counterclockwise 
directon.  Positive  force  outputs  act  in  the  directions  indicated  in 
Figure  5.3.  The  accelerometers  to  be  used  as  sensors  were  mounted 
colocated  with  the  actuators  with  the  same  numbering  and  output  sign 
(see  Figure  5.3).  Accelerometer  measurements  were  input  to  the  signal 
conditioning  and  integration  circuit  to  estimate  beam  velocity.  This 
velocity  measurement  was  then  filtered  using  an  A/C  coupled  amplifier  to 
remove  any  bias  in  the  integration  circuit.  Finally,  the  measurements 
were  input  to  the  PC-1000,  which  is  a  high  speed  digital  array  processor 
used  for  real-time  data  collection,  estimation,  and  control 
applications.  Analog  inputs  are  converted  to  12-bit  digital  data  and 
internal  calculations  are  performed  using  32-bit  floating  point 
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arithmetic.  Calculated  outputs  are  converted  back  to  analog  signals  to 
command  the  actuators.  The  processor  can  be  programmed  for  various 
sample  rates  (in  Hz),  input  and  output  signal  ranges  (in  volts), and 
input/output  channel  gains.  Further  information  about  the  PC- 1000,  its 
operation  and  capabilities,  can  be  found  in  Reference  12.  The  PC- 1000 
was  programmed  with  the  estimation  and  control  algorithm  to  be  tested, 
and  the  control  commands  it  generated  were  input  to  the  actuators.  In 
operating  the  PC-1000,  a  small  bias  was  found  on  the  input  channels, 
probably  residing  in  the  analog- to-digital  (A/D)  converter.  This  bias 
would  pass  through  the  estimator/controller  and  saturate  the  actuators 
To  remove  it,  the  control  commands  were  filtered  using  a  set  of  A/C 
coupled  amplifiers  identical  to  those  used  in  filtering  the  sensor 
measurements.  The  conditioned  control  commands  were  then  provided  to 
the  actuators. 


Svs tern  Gain  Verification 


Prior  to  attempting  to  control  the  structure  with  the  PC- 1000,  an 
analog  velocity  feedback  loop  was  implemented  to  verify  the  gains 
through  the  system.  Using  the  sensor  scale  factors  calculated  in 
Chapter  4  and  the  actuator  transfer  functions  from  Chapter  2,  Table  2.7, 
a  viscous  damping  coefficient  in  Ibf/in/sec  can  be  calculated.  The 
coefficient  can  be  changed  using  the  Intek  amplifiers  in  the  sensor 
signal  conditioning  path,  which  could  be  adjusted  in  10  dB  increments. 

Before  closing  the  feedback  loop,  the  open  loop  damping  due  to  the 
actuators  was  characterized  by  finding  a  damping  coefficient  which 
would  generate  a  structure  frequency  response  that  matched  the  damping 
measured  in  the  final  modal  test.  This  was  done  using  the  same  method 
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discussed  in  Chapter  2  for  estimating  the  acutator  residual  damping  (see 
equation  5  -  7) .  This  residual  damping  coefficient  was  then  added  to 
the  commanded  feedback  coefficient  for  each  actuator  controlling  the 
axis  under  test  resulting  in  a  total  coefficient  applied  to  the 
structure.  Two  tests  were  conducted  in  each  bending  axis,  one  with  the 
nominal  sensor  gains  calculated  in  Chapter  4,  and  a  second  with  the 
sensor  scale  factor  increased  by  10  dB.  Table  5.1  summarizes  the 
damping  coefficients  calculated  for  the  closed  loop  response  for  these 
tests. 


Table  5.1.  Closed  Loop  Viscous  Damping  Coefficients  (Ibf/in/sec) 


These  damping  coefficients  were  used  in  the  finite  element  model  to 
generate  a  frequency  response  for  the  structure.  The  damping  ratio 
predicted  for  the  structure  was  estimated  from  the  frequency  response 
using  the  half  power  method  (equation  8).  Tests  were  conducted  on  the 
structure  by  displacing  the  base  plate  and  recording  the  free  decay  time 
response  of  the  beam  on  the  Ono  Sokki  FFT  Analyzer.  Since  the  Ono  Sokki 


records  the  data  digitally,  the  peaks  of  the  response  could  be  listed 
and  the  damping  ratio  estimated  using  the  log  decrement  method  (7:61) 


where  M 


Magnitude  of  the  first  peak 


■  Magnitude  of  the  next  peak  in  the  response 


To  obtain  an  accurate  estimate,  f  was  calculated  for  each  pair  of  peaks 
in  the  response.  The  damping  ratio  and  uncertainty  is  listed  for  each 
test  in  terms  of  the  mean  and  standard  deviation  of  all  damping  ratios 
calculated  for  a  single  free  decay  of  the  structure. 

Table  5.2  summarizes  the  predicted  and  measured  damping  ratios  for 
the  experiment.  Figures  5.4  and  5.5  show  the  predicted  frequency 
response  and  the  measured  time  response  for  the  z  axis  tests,  and 
Figures  5.6  and  5.7  are  the  same  data  for  the  y  axis  tests.  Agreement 
between  predicted  and  measured  damping  ratios  is  excellent  for  these 
tests.  The  slight  differences  are  most  likely  due  to  the  damping  added 
by  the  sensor  and  actuator  wiring,  not  accounted  for  in  the  model. 


Table  5.2.  Predicted  vs  Measured  Damping  Ratios 


Test  » 


Predicted  f 


0.0590 


0.0722 


Measured  f 


0.059 


0.0424 


0.0579 


0.074 


0.044 


0.061 
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Figure  5.5,  Predicted  and  Measured  Z-Axis  Response  (C  -  0.0565) 


With  the  results  of  the  analog  closed  loop  tests  indicating  the 
sensor  and  actuator  scale  factors  are  correct,  the  PC-1000  input  gains 
were  calculated  using  the  accelerometer  scale  factors.  These  gains  are 
programmed  into  the  PC- 1000  and  boost  the  sensor  scale  factors  to  unity 
prior  to  the  signal  going  to  the  estimator/controller .  Gains  for  each 
channel  are  listed  in  Table  5.3. 

Table  5.3.  PC-1000  Sensor  Input  Gains 


Sensor  Channel  I  Gain  (Ibf/ft/sec)  Gain  (Ibf/in/sec) 


1 


10.847 

130.208 

10.723 

128.700 

11.099 

133.156 

11.099 

133.156 

VI .  Closed  Loop  Digital  Control 


Theory  and  Implementation 


In  general,  the  equations  of  motions  for  a  vibrating  structure 
which  is  being  controlled  can  be  written  in  the  form 


Mq  +  Cq  +  Kg 


where  g  is  an  n-vector  of  generalized  coordinates,  M  is  an  nxn  symmetric 
mass  matrix,  C  is  and  nxn  damping  matrix.  K  is  an  nxn  stiffness  matrix, 
u  is  an  m-vector  of  control  inputs  to  the  structure,  and  D  is  a  nxm 
matrix  of  actuator  coefficients  (Ref  4).  Equation  16  can  be  decoupled 
',:£ing  the  modal  coordinates  r)  aiid  the  modal  matrix  of  right  eigenvector.s 
,  to  define  a  transformation  such  that 


wh.ich,  after  substituting  for  g  and  its  higher  time  derivatives,  and 
7  remul  t  i p  1  V ing  bv  <t>‘,  transforms  equation  16  to  the  form 

"  *2  T 

r)+[2fu)]r7+[u)]i7  ■=  'I’  Du  ( 

wliere  <t>  has  been  normalized  to  meet  the  criteria  =  I,  an  nxn 


identity  matrix,  !2ru']  is  an  nxn  diagonal  damping  matrix  and 


w  IS  an 


nxn  diagonal  matrix  of  the  eigenvalues  of  equation  16,  Equation  18  can 
be  formulated  in  the  standard  state  space  form  of 


6  -  1 


‘ j  vjv  *->  , ' .  -  > .  1  k .  -v -N -•Si  V. 


(19) 


X  -  Ax  +  Bu 

where  x  is  defined  as  a  2n-vector  composed  of  r)  and  and  A  and  B  have 
the  form 
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(20) 


In  general,  the  states  of  a  structure  cannot  be  directly  measured, 
but  are  some  linear  combination  of  the  generalized  coordinates  g  which 
can  be  measured,  thus  equation  19  is  supplemented  with  the  measurement 
equat ion 

y  -  C  3  ( 21 ) 

where  v  is  an  s-vector,  s  being  the  number  of  sensors  used  in  the 
system,  and  C  is  an  sxn  matrix  which  can  be  partitioned  into 

C  -  [  C  I  C  ]  (22) 

where  the  partitions  C  and  are  coefficients  for  position  and 
velocity  sensors  respectively.  Using  the  modal  matrix  to  write  equation 
21  in  terms  of  the  states  gives  the  form 


1 


[  C  <!>  1  C  ] 


X 


(23) 


The  state  space  model  for  the  system  is  represented  by  equations 
19  and  23.  For  this  experimental  configuration,  C  will  be  a  zero 

p 

matrix  since  only  velocity  is  being  derived  from  the  accelerometer 
measurements.  Closing  the  control  loop  using  state  space  techniques 
simply  requires  feeding  back  the  states  through  a  gain  matrix  G 
(11:327).  Thus  the  control  input  is  defined  by 

u  -  -G  X  ( 2-^  ' 

Unfortunately,  since  the  controlled  states  cannot  be  measured  directlv. 
an  estimator  is  required  to  generate  the  states  to  be  used  with  the 
control  matrix.  One  method  of  constructing  this  observer  is  to  use  the 
system  model  to  estimate  the  states  based  on  the  sensor  measurements. 

For  this  formulation,  the  estimator  has  the  form 

• 

X  “  Ax  +  Bu  +  K(y  -  y) 
y  -  Cx 

whore  x  is  the  estimated  state  and  y  is  the  estimated  output.  The 
observer  gain  matrix  K  is  chosen  such  that  the  error  in  the  state 
estimate,  defined  by 

e  -  X  -  X  (27 ) 

is  stable.  Since  control  must  be  based  on  the  estimated  states,  the 
control  input  u  must  have  the  form 

u  -  -G  X  (28) 


(25) 

(26) 
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Substituting  for  y  and  u  in  equation  25  and  collecting  like  terms 


results  in 

X  -  (A  -  BG  -  KC)x  +  Ky  (29 

which  with  equation  28  defines  the  estimator/controller  to  be  used  in 
the  feedback  loop. 

There  are  many  methods  for  selecting  the  K  and  G  matrices  to  give 
the  plant  the  desired  characteristics.  The  one  that  will  be  used  for 
this  experiment  is  the  optimal  linear  quadratic  regulator.  This  method 
determ.ir.es  G  bv  minimizing  the  performance  index 


w}iere  Q  and  R  are  weighting  matrices  chosen  such  that  the  eigenvalue.s  o 
('A  -  BG)  ,  the  controlled  plant,  exhibit  the  desired  stablity.  K  is 
similarly  determined  by  minimizing  the  performance  index 
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v*R  v)  dt 
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wherc  Q^^and  R  ^  are  again  weighting  matrices  selected  such  that  the' 
eigenvalues  of  (A  -  KC j  define  an  observer  with  the  desired  stablitv. 
For  this  experiment,  R  and  R  ^  were  selected  as  identity  matricc.s  aixi  Q 
and  Q  ^  were  selected  independent  1  v  to  give  the  desirt-d  estimator  and 
controller  eigenvalues. 

To  implement  the  estimator/controller  on  the  Pd-ltUili,  tht  solvit  iot 
to  equation  25  is  discretized  into  time  steps  Defining 


the  solution  for  a  discrete  time  step  6t  becomes  (5;25) 


exp(A  fit)  X  +  A  ^[exp(A  fit)  -  I]K  y 

oc  k  oc  oc  “•  •‘'k 

(33) 

where  fit  is  defined  by  the  inverse  of  the  sample  rate  used  in  the  array 
processor.  The  computational  format  of  the  PC-1000  is  (Ref  12) 
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where  the  dimensions  of  the  partitions  are  F  ;  16x16,  F  : 16x32, 

^  11  12 

F  :32xl6,  and  F  :32x32.  Puttine  equations  33  into  this  form  defines 
2;  22  0-1 

the  partitions  of  the  F  matrix  to  be 


F  -  0 

11 

F*-G 

12  (35) 

F  -  A'Hexp(A  fit)  -  iIk 

21  oc  ^  oc  ^ 

F  -  €xp(A  St) 

22  o  c 

The  specific  A,  B,  and  C,  matrices  used  to  define  the  plant  are 
listed  in  Appendix  F.  Calculation  of  the  estimator/controller  gain 
matrices  and  the  system  eigenvalues  was  accomplished  using  software 
based  on  work  previously  accomplished  by  Aldridge  (Ref  13)  and  operating 
on  one  of  the  Vax  11/785  computers  resident  at  AFIT.  The  gain  and 


system  matrices  were  input  to  Matlab,  a  matrix  manipulation  software 
package  residing  on  the  Vax .  which  was  used  to  calculate  the  F^,  and 
matrices . 


Closed  Loop  Simulation 


A  closed  loop  simulation  was  implemented  on  the  PC- 1000  to  verify 
the  time  response  of  the  closed  loop  system  is  an  accurate  indicator  of 
the  predicted  closed  loop  plant  eigenvalues,  and  to  evaluate  the 
processor's  ability  to  perform  the  estimator/controller  function.  The 
estimator/controller  configurations  formed  for  the  simulation 
incorporated  the  first  three  natural  modes  for  the  structure.  In 
formulating  the  controller  gain  matrix  the  second  and  third  modes  were 
deweighted  in  the  cost  function,  thus  control  was  applied  only  to  the 
first  2-axi.s  bending  mode.  Observer  gains  were  weighted  approximately 
equally.  Two  controllers  were  generated  with  identical  control  gains 
and  different  observer  stablities  to  predict  the  effect  of  estimator 
errors  on  the  closed  loop  system.  Slave  processor  =2  of  the  array 
processor  was  loaded  with  the  open  loop  plant  model  (x  =  Ax  -f  Bu)  while 
slave  processor  ~1  remained  as  the  estimator/controller.  An  initial 
condition  of  x^  -  0.01  was  placed  in  the  plant  model  which  approximates 
an  initial  displacement  of  the  base  plate  in  the  z  direction.  The  time 
response  of  various  parameters  in  the  closed  loop  system  was  used  to 
estimate  the  actual  damping  expected  in  the  first  z  axis  bending  mode 
when  control  is  applied  to  the  beam.  Data  was  recorded  on  the  Ono  Sokki 
spectrum  analyzer  and  damping  was  estimated  using  the  log  decrement 
method  (equation  15)  on  the  time  response  data.  Comparisons  of  the 
actual  and  estimated  states  and  the  sensor  output  and  actuator  commands 


are  shown  in  Figures  6.1  through  6.6.  A  summary  of  the  system  damping 
estimated  from  the  various  parameters  in  the  system  is  presented  in 
Table  6.1.  Estimator  damping  and  modal  damping  values  are  obtained  from 
the  eigenvalues  of  (A  -  KC)  and  (A  -  BG)  respectively. 


Table  6.1.  Simulation  Damping  Results  for  First  Z  Bending  .Mode 


Estimator 

Modal 

Plant 

States 

Sensor 

Estimated  States 

Damping 

Damping 

^  X 

1 

^  s 

2 

1 

. 

0.155 

0.101 

0.093 

0.111 

0.113 

0.068 

0.091 

0.548 

0.101 

0.091 

0.101 

0.101 

0.081 

0.099 

Comparing  the  plant  damping  with  the  predicted  modal  damping 
indicates  tiie  less  damped  estimator  should  overcontrol  the  plant 
slightly.  This  is  probably  due  to  the  error  dynamics  in  the  estimated 
states  increasing  the  modal  amplitude.  As  estimator  damping  increases, 
the  error  has  less  time  to  contribute  to  the  estimated  state  and  the 
plant  damping  approaches  the  predicted  modal  damping.  This  is  readily 
seen  in  Figures  6.7  and  6.8  which  show  the  estimator  response  to  an 
initial  condition  with  the  sensor  input  removed  from  the  PC-lOOf'. 
Another  observation  is  that  the  sensor  data  is  a  more  accurate 
measurement  of  the  actual  plant  damping  than  the  estimated  states,  even 
for  a  highly  damped  estimator. 
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Command  Output  -  Estimator  C  -  0.5^8 


SEC 


Command 


Closed  loop  testing  used  the  integrated  accelerometer  measurements 


from  the  beam  as  velocity  inputs  to  the  estimator  and  the  control 
commands  from  the  PC-1000  were  input  to  the  proof  mass  actuators.  The 
estimator  was  started  with  zero  initial  conditions  and  the  beam 
quiescent.  Even  so,  transients  proved  to  be  a  problem.  The  PC- 1000  was 
found  to  have  a  small  bias  voltage  on  the  input  channels  (approximatelv 
60  mV),  probably  in  the  A/D  converter.  This  bias  was  sufficient  to 
cause  high  gain  controllers  to  produce  transients  large  enough  to 
saturate  the  actuator  stroke.  To  prevent  overdriving  and  possiblv 
damaging  the  actuators,  power  was  removed  until  all  transier“s  settled. 

Controllers  were  designed  for  the  fundamental  bending  mode  in  each 
direction,  and  the  first  torsion  mode.  Estimator/controller  states 
included  all  three  modes,  however  control  weightings  were  used  which 
applied  control  to  only  one  mode  at  a  time  and  left  the  remaining  modes 
essentially  unchanged.  Structural  response  for  the  fundamental  modes 
was  measured  by  giving  the  beam  an  initial  displacement  or  velocity  in 
the  controlled  axis  and  measuring  the  time  response  of  sensors  on  that 
axis.  Damping  ratios  were  calculated  for  each  test  using  the  log 
decrement  method.  Results  from  several  estimator/controller 
configurations  are  showm  in  Table  6.2. 

The  measured  damping  coefficients  in  Table  6.2  indicate  a  trend 
opposite  to  that  seen  in  the  simulation.  As  the  ratio  between  thf 
estimator  and  controller  decreases,  the  measured  damping  also  decreases 
The  simulation  predicted  higher  damping  in  the  case  where  the  estimator 
damping  closely  approached  the  modal  damping. 

To  determine  if  the  log  decrement  data  was  accurate,  the  GenRad 
modal  analyzer  was  also  used  to  measure  the  modal  damping  for  the  first 
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Table  6.2.  Closed  Loop  Modal  Damping  Coefficients 


Estimator 

Damping 

Predicted 

Damping 

Ratio 

Measured 

Damping 

0.3576 

0.0716 

4.994 

0.080  ±  0.020 

0.2798 

0.0716 

3.908 

0.069  ±  0.018 

0.2149 

0.0716 

3.001 

0.070  ±  0.009 

0.1214 

0.0716 

1.696 

0.066  ±  0.013 

0.1015 

0.0716 

1.418 

0.066  ±  0,020 

0.0793 

0.0716 

1.112 

0.041  ±  0.016 

Error 

(%) 


0.1523  0.0715  2.130  0.059  ±  0.011 
0.1252  0.0715  1.751  0.058  ±  0.015 
0.0980  0.0715  1.371  0.056  ±  0.007 
0.0865  0.032A  2.670  0.030  ±  0.010 


two  bending  modes.  The  beam  was  excited  in  the  controlled  axis  using  a 
force  hammer  and  damping  estimates  were  obtained  using  a  complex 
exponential  curve  fit  to  the  resulting  frequency  response  function 
measured  by  the  GenRad.  A  comparison  of  these  results  to  the  log 
decrement  data  is  shown  in  Table  6.3.  While  the  curve  fit  data  appears 
to  be  more  accurate  and  shows  the  trend  more  clearly,  the  damping 
estimate  was  highly  dependent  upon  the  frequency  range  the  software  was 
allowed  to  fit.  For  the  z-axis  test  with  the  estimator  at  0.155, 
damping  estimates  ranging  from  0.088  to  0.056  were  observed.  The  data 
presented  represents  a  "best  fit”  to  the  frequency  response  function 
based  on  minimizing  the  resi-^uals  between  the  measured  and  fit  data. 
Since  the  log  decrement  data  is  reasonably  accurate,  testing  of  the 
fundamental  bending  and  first  torsion  modes  continued  using  this  method 


Table  6.3.  Damping  Estimation  Method  Comparison  -  Log  Decrement  vs 

Complex  Exponential  Curve  Fit 


Axis 

Estimator 

Damping 

Modal  Damping 
(predicted) 

Modal  Damping 
(log  dec) 

Modal  Damping 
(curve  fit) 

z 

0.1214 

0.0716 

0.066 

0.077 

z 

0.1015 

0.0716 

0.066 

0.063 

Z  1 

0.0793 

0.0716 

0.041 

0.043 

y 

0.1252 

0.0715 

0.058 

0.059 

In  an  attempt  to  determine  a  cause  of  the  lower  than  expected  modal 
damping  several  more  tests  were  run  with  the  simulation  to  measure  gains 
through  the  controller/estimator.  Table  6.4  shows  the  gains  between 
variou.s  system  parameters. 


Table  6.4.  System  Gains  (dB) 


Estimator  ^ 

X  /x 
r  1 

X  /s 
r  2 

• 

X  /s 
r  2 

u  /s 

2  2 

0.155 

-3.0 

-2.8 

-34.9 

-16.5 

-29.2 

0 . 548 

-2.2 

o 

00 

-34.9 

-14.7 

-27.5 

Predicted  modal  damping  of  0.101  for  both  controllers 

The  simulation  shows  that  the  sensor  to  acutator  command  gain  was 
smaller  for  the  less  damped  estimator.  This  result  is  consistent  with 
tlie  trend  in  the  experimental  data  of  less  modal  damping  for  the  less 
damped  estimators,  and  is  also  consistent  with  the  fact  that  the  higher 
damped  estimators  have  larger  gains  in  the  estimator  gain  matrix  K.  The 
discrepancy  which  remains  unresolved  is  that  the  simulation  still 
maintains  the  predicted  level  of  modal  damping  even  though  it  is 


providing  less  control  input  to  the  plant.  Further  investigation  of  the 
simulation  would  be  required  to  fully  evaluate  this  difference. 

A  tabulation  of  additional  single  axis  control  tests  conducted  is 
in  Table  6.5.  From  this  data,  to  have  the  actual  modal  damping  be  close 
to  the  predicted  eigenvalue,  the  estimator  damping  should  be  at  least 
twice  the  modal  damping  for  z-axis  control,  and  three  times  the  modal 
damping  or  higher  for  y-axis  and  torsion  control. 

While  keeping  the  estimator  damping  high  relative  to  the  modal 
damping  is  desirable,  this  was  found  to  limit  the  modal  damping  level 
that  could  be  used  in  the  experiment.  The  estimators  operating  with  a 
damping  coefficient  greater  than  0.2  were  found  to  pass  low  frequency 
signals  through  the  controller.  When  the  control  gains  are  large,  the 
low  frequency  drift  in  the  sensor  signal  caused  the  actuators  to  wander. 
For  some  controllers  the  drift  was  so  bad  that  the  actuators  would  not 
remain  centered  long  enough  to  run  even  a  single  decay  test.  Figures 
6.9  through  6.12  show  the  actuator  command  versus  the  sensor  input  for 
different  estimator  damping  ratios.  A  z-axis  controller  with  an 
estimator  damping  coefficient  of  0.356  was  the  highest  that  could  be 
tested.  Estimators  operating  with  damping  coefficients  less  than  about 
0.2  had  very  little  problem  with  drift.  While  some  was  still  present, 
the  magnitude  was  either  below  the  actuator  response  threshold  or  slow 
enough  to  allow  adjustment  of  the  actuator  centering  before  a  test. 


Predicted 


Experimental 


Estimator 

Damping 

Modal 

Damping 

Ratio 

f'/r 

X  X 

0.1015 

0.0716 

1.9176 

0.1268 

0.0999 

1.2757 

0.1590 

0.1295 

1.2369 

0.1863 

0.1507 

1.2369 

0.2199 

0.1790 

1.2006 

0.1268 

0.0716 

1.7709 

0. 1590 

0.0999 

1.5993 

0.1863 

0.1295 

1.9969 

0.2199 

0.1507 

1.9260 

0.0980 

0.0715 

1.3706 

0.1252 

0.0965 

1.2979 

0.1523 

0.1253 

1.2155 

0.1866 

0.1992 

0.2136 


0.1252 


0.1523 


0.1866 


0.2136 

0.2^58 


0.1238 
0  .  U92 
0.1782 


1.5073 
1.9316 
1  .  3793 


Measured 
Modal  Damping 


%  of 
Predicted 


0.066 

0.020 

0.085 

+ 

0.013 

0.105 

+ 

0.022 

0.128 

+ 

0.031 

0.159 

4 

0.065 

0.065 

4 

0.025 

0  081 

4 

0.013 

0.115 

4 

0.030 

0.138 

4 

0.093 

0.056 

4 

0.007 

0.076 

4 

0.017 

0.099 

4 

0.019 

0.119 

4 

0.009 

0.193 

4 

0.009 

0.058 

4 

0.015 

0.081 

4 

0.015 

0. 102  1  0.010 
0.125  ±  0.006 
0. 153  +  0.035 


0.922  ±  0.279 


0.855  ±  0.131 


0.893  ±  0.177 


i 


0.815 

4 

0.131 

0.929 

4 

0.291 

0.916 

4 

0.285 

0.783 

4 

0.098 

0.786 

4 

0.176 

0.790 

4 

0.152 

0.769 

4 

0.027 

0.802 

4 

0.022 

0.811 

4 

0.208 

0.839 

4 

0.155 

0.829  ±  0.081 
0.838  1  0.090 
0.859  ±  0.196 


Output  Drift  -  Sensor  Input  vs  Command  Output  (C  -  0.3576) 


Sensor  Input  vs  Command  0i 


Sensor  Input  vs  Command  Output  ( (T  -  0.21^49) 


VII.  Modal  Su 


Theory 

Using  a  reduced  order  dynamic  model  reduces  the  order  of  the 
controller,  and  thus  the  computational  requirements  can  be  significantly 
reduced  in  a  very  large  system.  If  the  reduction  is  accomplished  bv 
just  deleting  modes,  the  uncontrolled  modes  can  be  driven  unstable. 

This  occurs  because  the  sensor  outputs  and  control  commands  still 
contain  information  about  the  deleted  modes.  These  effects  have  becomt 
commonly  refered  to  as  observation  and  control  spillover,  A  method  for 
eliminating  spillover  and  maintaining  a  stable  system  was  developed  by 
Calico  and  Janiszewski  (Ref  A),  and  is  used  in  this  experiment  to 
construct  the  reduced  order  controllers. 

To  accomplish  this,  the  modes  of  a  system,  are  classified  as 
controlled,  suppressed,  and  residual.  The  state  vector  for  the  system 
now  takes  the  form  of 


x-(x,x.x>'  (3('!' 

-  —c  ~s  — r 

where  x  is  an  n  -vector  of  controlled  states,  x  is  an  n  -vector  of 

—  c  c  — s  s 

suppressed  states,  and  x  is  an  n  -vector  of  residual  states.  The  moder. 
included  in  x  are  onlv  those  necessary  to  establish  satisfactory  system 

—  C  ' 

performance  and  do  not  necessarily  include  the  lowest  order  modes. 

Using  this  partitioning,  the  state  space  model  can  be  written  as 


X 

c 

-  A 

C 

X  -t  B  u 
—  c  c  — 

(37  . 

• 

X 

-  A 

X  +  B  u 

V  38;) 

—  s  s  ~s  s  — 
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(39) 


X  -  A  X  +  B  u 

—r  r  —r  r  — 

y-Cx+Cx+Cx  (90) 

^  c  — c  &  — s  r  — r 

The  control  input  is  still  defined  by  equation  24,  but  with  x  replaced 
hv  X  .  The  coupling  of  the  control  u  into  the  suppressed  and  residual 
equations  could  excite  and  destabilize  these  inodes,  and  as  such  the 
problem  is  known  as  control  spillover. 

A  similar  situation  exits  for  the  the  observer.  The  estimator 
still  has  the  form  of  equations  25  and  26, 


X  -Ax  +Bx  +K(v-y) 

C  C  C  C  c  •- 

(41) 

y  -  C  X 

c  c 

and  as  such,  v  still  includes  information  about  the  suppressed  and 
residual  states.  This  coupling  can  induce  errors  into  the  estimated 
states  and  thus  generate  inappropriate  control  commands,  possibly 
driving  the  system  unstable.  This  is  the  effect  known  a  observation 
sp i 1  lover . 

To  maintain  system  stability,  it  is  sufficient  to  eliminate  either 
tvpe  of  spillover  (Ref  4).  For  this  experiment  it  was  decided  to 
eliminate  observation  spillover,  which  can  be  accomplished  by 
constraining  the  estimator  gain  and  output  matrices  such  that 

KC  0  (42a) 

C 

KC  -  0  (42b) 

s 


For  a  reduced  order  controller  the  residual  modes  are  ignored  and  the 
solution  to  equation  42b  can  be  found  by  singular  value  decomnos i t i on  of 


C^.  If  is  of  full  rank,  a  solution  exists  only  if  the  number  of 
sensors  is  greater  than  the  number  of  suppressed  modes .  This  is 
necessary  because  eliminating  the  observation  requires  constructing  a 
new  measurement  set  w,  which  is  a  linear  combination  of  sensor  outputs 
such  that  the  suppressed  mode  has  in  effect  been  subtracted  out  of 
the  measurement,  thus  suppressed.  This  in  effect  "costs"  a  sensor  for 
a  mode  to  be  suppressed.  The  results  of  the  singular  value 
decomposition  of  C  are  an  orthogonal  set  of  left  singular  vectors  which 
can  be  partitioned  into  a  set  associated  with  the  non-zero  singluar 
values  of  C  ,  and  a  set  associated  with  the  zero  singular  values  of  C^. 
Defining  the  set  associated  with  the  zero  singular  values  as  f,  equation 
42b  becomes 


rc  -  0  (43) 

S 

Defining  a  new  relation  to  represent  the  output  after  the  suppressed 
modes  are  removed 


w  - 


(44) 


equation  40  becomes 


w  -  rc  X  -r  rc  X  +  rc  x 

—  c  — c  5  r  ~r 


which,  by  ignoring  residuals  reduces  to 


w  -  rc  X 

—  c  ~c 


(45) 


(46) 
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using  this  new  measurement  in  the  estimator  equation  results  in 


X  -  A  X  -  B  u  -  K(w  -  w) 

c  c  c  —  — 


w  -  rc  X 

c  c 


u  -  -Gx 


which,  when  put  in  terms  of  defines  a  state  estimator  of  the  form 


X  -  (A  -  B  G  -  KrC  )x  +  KFv 

c  c  c  c  C  C 


The  time  domain  solution  of  equation  48  and  its  formulation  for  use  in 
the  PC- 1000  is  identical  to  equations  32  through  35,  when  the  full  state 
matrices  are  replaced  with  their  reduced  order  counterparts,  and  the 
estimator  gain  matrix  K  is  replaced  with  KF . 


Results 


Since  the  first  two  modes  of  the  system  are  the  first  bending  modes 
in  two  orthogonal  directions,  they  are  already  physically  decoupled  and 
it  is  not  necessary  to  apply  modal  suppression  before  eliminating  the 
orthogonal  axis  from  the  controller.  Therefore,  suppressed  mode 
controllers  were  designed  and  implemented  for  the  z  axis  of  the  beam 
which  dealt  with  the  first  and  second  z  bending  modes  and  torsion 
Initially,  to  determine  the  effect  of  spillover  on  the  system,  a  two 
mode  controller  using  only  torsion  and  second  z  bending  was  constructed, 
with  the  first  z  bending  included  in  the  model  as  a  residual  mode  to 
predict  its  stability.  The  predicted  eigenvalues  for  the  estimator  and 


I 
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controller  without  incorporating  the  residual  mode  are  listed  in 
Table  7.1.  The  eigenvalues  for  the  system  with  the  residual  mode 
incorporated  are  shown  in  Table  7.2. 


Table  7 . 1 . Predicted  System  Eigenvalues  Without  Residual  Modes 


Mode 

Estimator 

Controller 

r 

U> 

u> 

T1 

0.1029 

79.8805 

0.0126 

79.0245 

Z2 

0.1055 

11A.974 

0.0301 

113.688 

Table  7.2.  Predicted  System  Eigenvalues  With  Residual  Modes 


Mode 

— 

Estimator 

... 

Controller 

Residual 

— 

r 

u> 

a> 

T1 

0.1029 

79.8805 

0.0126 

79.0245 

Z2 

0.4701 

65.6812 

0.0048 

114.387 

Z1 

-1.0 

-1.0 

21.0/09 

10,0318 

While  the  system  appears  well  behaved  when  the  residual  mode  is 
ignored,  the  overall  system  is  obviously  unstable  when  the  omitted  mode 
is  accounted  for.  To  verify  this  prediction,  the  estimator/controller 
was  implemented  on  the  PC-1000.  The  resulting  control  response  would 
immediately  saturate  the  actuators  when  initiated,  and  eventually  drove 
the  z-axis  unstable. 
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To  overcome  the  loss  of  observability,  the  sensors  on  the  y  axis 
were  changed  to  stations  9  and  10  on  the  z  axis.  Modal  amplitudes 
listed  in  Table  3.7  indicate  the  response  at  these  stations  is  the  same 
order  of  magnitude  as  the  first  bending  mode  at  the  base  plate 
Therefore,  the  second  z  bending  signal  should  be  substantial  enough  to 
not  be  subtracted  out  when  the  first  mode  is  suppressed 

The  new  sensor  stations  were  already  instrumented  with  PCB 
Structural  Accelerometers  used  in  the  modal  tests.  Rather  than 
re  instrument ,  the  PCB  accelerometers  were  used  in  place  of  the  Kistlers 
The  new  accelerometer  calibrations  and  scale  factors  after  signal 
conditioning  are  shown  is  Table  7.5. 


Table  ,  5  PCB  Accelerometer  Calibrations  and  Scale  Factors 
Station  Serial  «  Calibration  (V/g)  Scale  Factor  (mV/ft/sec) 


93  .  10 


1223 

0.8902 

1216 

0.9321 

The  results  of  the  new  sensor  configuration  are  shown  in  Tables  7.6 
and  7.".  The  observabi 1 i tv  of  the  second  mode  returned  and  the 
.supf  ress  ion  of  the  first  bending  had  minimal  effect  on  the  controlled 
s vs  tern  eigenvalues. 


Table  7.6.  Predicted  Svstem  Eigenvalues  Before  Suppression  -  Adjusted 

Sensor  Locations 


Estimator 


Control  1 p  r 


Mode 

- 

U) 

r 

T1 

0.0876 

79.6438 

0.0127 

79.0425 

Z2 

0.0827 

114.473 

0.0151 

113.688 

Table  7.7.  Predicted  System  Eigenvalues  After  Suppression  -  Adjusted 

Sensor  Locations 


Mode 

Estimator 

Controller 

u 

C 

U' 

T1 

0.0876 

79.8412 

0.0126 

79.0425 

Z2 

0.0824 

114.467 

0.0151 

113.688 

Implementation  of  this  controller  proved  to  be  only  partially 
.successful.  The  first  z  bending  mode  was  not  driven  unstable,  however 
is  was  slightly  destablized.  Figure  7.1  shows  the  time  response  for  the 
first  z  bending.  The  damping  ratio  calculated  for  the  test  was  0.0267. 
a  red'jction  of  0.01  over  the  normal  residual  damping  in  the  structure. 
Ttrv  damping  on  the  second  z  bending  was  estimated  using  the  impact 
mi ti.od  and  the  modal  analyzer.  The  transfer  function  was  then  curve  fit 
to  determine  the  dan. ping,  the  result  being  shown  in  Figure  7.2.  As  is 
readily  apparent  from  the  estimated  roots,  the  damping  has  only  been 
marginally  increased  from  its  residual  value  of  0.00291  to  0.00993. 

'/iti,  a  predicted  controlled  damping  of  0.0151,  this  is  by  far  the 
largest  deviation  from  prediction.  The  most  likely  explanation  is  the 
low  modal  amplitude  of  the  second  mode  at  the  base  plate.  The  sen.sor 
input  and  actuator  command  for  this  controller  are  shown  in  Figure  7.3 

A  problem:  which  ran  be  seen  in  Figure  7  3  was  a  low  frequencv 
tr.TTiSient  in  the  PCB  response  when  the  beam  was  excited.  The  controllc)' 
output  followed  this  signal  verv  closely  and  almost  invariably  saturated 
the  actuator  stroke  unless  the  force  input  was  kept  extremely  small. 

Finally,  to  determine  how  well  the  suppression  of  the  first  bending 
mode  is  operating  the  frequency  spectra  of  the  sensor  input  and  the 
suppressed  measurement  w  were  measured  while  exciting  the  suppressed 
mode  (see  Figure  7.4),  The  sensor  input  magnitude  of  -44.8  dBV  should 


Figure  7.1,  First  Z  Bonding  Mode  Response  -  Suppressed 


FiRi'.re  7.'*.  Sensor  vs  Suppressed  Measurement  Spectral  Response 


VIII ■ Conclusions  and  Recommendations 


Conclusions 

Several  issues  important  to  future  design  and  testing  in  the  area 
of  large  space  structure  control  have  been  raised  in  this  experiment. 
First,  if  linear  proof  mass  actuators  are  to  be  used  in  controlling 
lightweight,  flexible  structures,  care  should  be  taken  to  ensure  that 
not  only  are  they  not  placed  at  or  near  nodes  on  the  structure,  but  that 
they  don't  create  a  node  at  the  control  point  by  their  mere  presence. 
Colocation  of  these  actuators  should  be  avoided  to  the  extent  possible. 
The  second  problem  associated  with  the  actuators  is  the  placement  of 
their  dynamics.  To  reduce  the  computational  burden  on  the  control 
system,  it  is  always  desirable  to  have  a  flat  response  from  the  actuator 
over  the  control  bandwidth.  To  accomplish  this,  along  with  the 
desirable  characteristics  of  maintaining  a  free  floating  and  centered 
proof  mass,  dictates  moving  the  low  order  poles  of  the  system  to  very 
low  frequency.  This  has  the  undesirable  effect  of  making  the  actuator 
very  responsive  to  drift  and  bias  shift.  Using  the  actuators  in  this 
configuration  requires  highly  accurate  instrumentation  and  excellent 
signal  conditioning  to  filter  out  all  signals  outside  the  control 
bandwidth . 

This  sensitivity  obviously  makes  sensors  an  important  issue. 
Integrating  accelerometer  data  to  generate  a  velocity  measurement  will 
continually  generate  problems  with  low  frequency  signals  unless  an 
accelerometer  capable  of  accurately  sensing  acceleration  frequencies 
below  1  Hz  and  amplitudes  to  less  than  a  tenth  percent  of  scale  is  used. 
None  of  the  piezoelectric  accelerometers  evaluated  for  use  in  this 


experiment  were  truely  adequate  for  the  task.  While  the  fundamental 
frequency  in  this  experiment  was  still  above  1  Hz,  a  large  structure  on 
orbit  may  have  many  modes  at  or  below  this  frequency.  The  only 
accelerometers  capable  of  adequately  performing  in  this  regime  may  be 
the  PIGA  (Pendulous  Integrating  Gyroscopic  Accelerometer)  type  widely 
used  in  inertial  navigation  platforms  or  others  of  similar  capability. 
Position  sensors  are  also  possibilities  for  use  in  the  control  design. 
While  no  position  sensors  were  used  or  evaluated  for  this  experiment, 
development  of  an  accurate  position  sensor  for  use  on  a  large  space 
structure  control  would  expand  the  options  available  to  the  control 
designer.  Benefits  of  a  particular  sensor  or  sensor  type  would  be 
dependent  on  the  specific  experiment  or  vehicle  to  be  controlled. 

With  the  number  of  problems  seen  with  noise  in  this  experiment, 
both  high  and  low  frequency,  a  deterministic  observer  may  not  be  the 
best  algorithm  with  which  to  approach  the  problem.  While  some  of  the 
difficulties  could  be  reduced  by  better  signal  conditioning  and  more 
accurate  instrumentation,  the  less  costly  solution  may  be  to  implement 
some  form  of  a  stochastic  estimator  in  the  controller. 

Lastly,  while  the  demonstration  of  reduced  order  controllers  was 
only  partially  successful,  all  indications  from  the  experiment  were  that 
it  is  a  viable  technique  given  adequate  observability  and  controlabi 1 i ty 
of  the  modes  in  question. 

Recommendations  for  Future  Work 

A  number  of  improvements  could  be  made  to  increase  the  accuracy  of 
the  experiment  and  perhaps  make  it  a  better  test  bed  for  evaluating 
control  methodologies.  Several  improvements  in  the  area  of  the 
actuators  could  be  readily  made.  The  actuator  shafts  should  be  replaced 


with  ones  which  have  a  much  harder  surface.  Scoring  of  the  current 
shafts  from  the  proof  mass  bearings  was  becoming  significant  even  before 
closed  loop  testing  of  the  beam  was  started.  Additionally,  the  lag 
prefilter  should  be  removed  from  the  compensator  circuit  and  the 
actuator  response  verified.  Model  predictions  indicate  the  phase  and 
magnitude  deviations  of  the  actuator  transfer  function  at  the  low  end  of 
the  control  bandwidth  were  due  to  this  prefilter  and  not  the  actuator 
dynamics.  As  a  third  step,  a  thorough  review  of  the  compensator  circuit 
for  ways  to  increase  the  accuracy  and  reduce  the  noise  in  the  system 
would  be  profitable.  Initially,  isolation  of  the  power  inputs  to  the 
operational  amplifiers  may  be  helpful. 

In  the  area  of  sensing,  substantially  better  low  frequency  signal 
conditioning  will  be  required  to  be  able  to  operate  a  deterministic 
estimator  with  a  damping  ratio  higher  than  about  0.25. 

The  bias  observed  in  the  input  channels  of  the  PC- 1000  should  be 
isolated  and  removed.  While  the  problem  can  be  worked  around,  it 
degrades  the  performance  and  dynamic  range  of  the  experiment.  For  some 
of  the  higher  gain  controllers  applied  to  the  higher  modes  the  output 
bias  was  observed  to  be  as  much  as  7  volts.  With  a  maximum  output  of  10 
volts,  the  dynamic  range  of  the  control  signal  has  been  significantly 
reduced . 

Enhancing  the  controlabil ity  of  the  higher  bending  modes  is  needed 
to  adequately  demonstrate  modal  suppression,  but  this  may  be  difficult 
to  accomplish  for  this  structure.  If  the  actuators  controlling  one  axis 
are  removed,  the  finite  element  model  predicts  only  a  few  percent 
increase  in  the  relative  amplitude  between  the  first  and  second  bending 
modes.  A  more  viable  approach  to  being  able  to  demonstrate  suppression 
of  coupled  modes  may  be  to  place  an  asymmetric  mass  on  the  base  plate 
large  enough  to  significantly  couple  the  bending  fundamental  modes. 
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Appendix  A:  Actuator  Component  Specifications  and  Configuration 
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Figure  A.l.  Linear  Motor  Specifications 
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HPD  SERIES-HERMETICALLY  SEALED  (PIN  TERMINATION) 


•  hermetically  sealed  BY  TIG  AND 

EB  WELDING 

•  IMPERVIOUS  TO  hostile  ENVIRONMENTS 

•  THROUGH.BORE  CONSTRUCTION 

HPD  Series  unttsare  similar  to  the  DC-D  and  HCD  Senes. 
Tunpiien  men  ^as  (TiQ)  ana  electror^  beam  (EBi  «veio 
rng  pfonoe  hemnetic  sealing  rnat  ii  tree  from  Qxidatron- 


proOupng  faults  that  may  cause  leakage.  For  this  reason. 
HPD  Senes  LVDT's  are  impervious  to  dm.  water,  steam 
spray,  and  most  corrosives  They  have  been  Qualified  at 
preuures  UP  to  1000  psig  C70  barsl  and  are  suitable  for 
numerous  high.preuure  applications.  HPD  units  employ 
a  giasS'Scaied.  pm-terminal  header  that  allows  the  core 
ano  core  rod  to  oass  through  the  unit.  HPD  units  have 
oouoie  magnetic  shielding  that  rrsakes  them  insensitive 
to  CMternsi  magnetic  influences. 


GENERAL  SPECIFICATIONS 

moot  .  ..i15V  DC  inommaO.  i  20  mA 


moot  ....  i15  V  DC  inommaO. 

Ooeratmg  Temper¬ 
ature  Range  .  .  0*F  to  *r160*F(-ie*Cto-»70*C) 

Survive  Temper 

iiure  Range  .  .  -6S*F  to  '•200*F 
(-55*C  io*95*Cl 
Null  VoJiage  .  .  .  0  V  DC 

RiPO'e  .  .  .  Less  tnjr  25  rnV  rms 

Lines'  tv  2  0.25%  fu'i  rang* 

StaO'i^tv  0.125%  full  scale 


Temperature 
Coefficient  of 
Scale  Factor  .  . 
ShocK^  Survival  .  .  . 
Vroratfon  Tolerance 
Coil  Form  Material 
Housing  Material 
Electrical 

Termination  .  . 
Output  impeoance 


0.04%/*F  (D.08%'*C) 

250  p  f or  11  milliseconds 
IDfuP  to  2  kHz 
High  Density,  glass-filled  poiynser 
AlSl  400  sanes  stainless  steel 

S-pm  terminal  heaaer 
Less  than  i  Onm 


PERFORMANCE  SPECIFICATIONS  AND  DIMENSIONS 


IVO’ 

aOMiKAi 

SCi 

M0D£. 

iiiiiiAa 

sa: 

aUMIIA 

iAACt 

•M  MPC 

:0  8sc 

lUMC 

so  124 

10 

mnpc 

:C7« 

4C 

^  Mi  Heo 

m 

n 

1MI  Nec 

tl  ooc 

u 

teaiMc 

S2DK 

4C 

MM  MO 

Siaoc 

IJ 

MM  MO 

20 

immmo 

S10JO 

t  0 

euroNst 


ORDERING  INFORMATION 

(fota  oj\  page  32  for  mitfucuoni  on  non  ro  use  ma  cn#rr  y 


DMEHSIONS 
I  iC*rV' 


awt 

Cm 

IMAM 

iMm 

iKftn 

u 

7 

2  40 

0  M 

0  55 

«4 

) 

3.23 

1  ifi 

DM 

4T 

4 

4  ID 

1  ac 

1-M 

77 

1 

471 

700 

2J3 

115 

10 

104 

310 

3  32 

1fi$ 

13 

1142 

4JC 

4  04 

tr 

14 

Ut2 

t.7t 

7» 

7M 

17 

7014 

120 

IK 

MC 

24 

74  47 

12  00 

U  SI 

weuT 

OUTfUT 

HOC 

Moc  «•: 

I  Mac 

e««*  '  tm  wfi 


DU 

JOS 


.livot _  ^  * 

r  I 

COMMON  '0 


•  at  UHC  tl  tSiaWvi* 

m;  X  esa**  rMnf« 


Figure  A  6.  LVDT  Spec  i  f  icat  ion.s 


1)  Pure  integrators  in  the  either  the  forward  or  feedback  paths 
exacerbated  the  low  frequency  problems  to  the  point  that  the 
configuration  was  impossible  to  test.  The  proof  mass  continually 
wandered  into  the  end  brackets  and  would  not  maintian  a  stable 
oscillation.  This  effectively  eliminated  the  option  of  implementing  an 
inverse  plant  model  to  cancel  the  low  frequency  dynamics  of  the  system. 


2)  High  frequency  noise  was  a  problem  thoughout  the  development. 
The  high  frequency  dynamics  in  the  open  loop  response  proved  to  be  a 
problem  when  the  acceleration  feedback  loop  was  added,  especially  when 
the  gain  was  turned  up  in  that  loop.  High  gain  in  the  acceleration 
feedback  loop  is  desireable  since  it  has  the  two  benefits  of  moving  the 
actuator  dynamics  lower  and  increases  the  ability  of  the  proof  mass  to 
float  freely  when  uncommanded.  When  the  gain  was  turned  up  to  a  high 
enough  level  to  effect  the  system  dynamics,  high  frequency  noise  would 
feed  back  into  the  motor  and  drive  the  system  unstable.  The  only 
solution  found  to  this  problem  was  to  low-pass  filter  the  acceleration 
signal  in  an  attempt  to  reduce  the  magnitude  of  this  noise.  This  was 
successful  to  some  degree,  but  adding  enough  filtering  at  a  frequency 
low  enough  to  remove  the  150  hz  dynamics  destroyed  the  phase  of  the 
feedback  signal  in  the  operating  bandwidth.  Three  low-pass  filters  with 
break  points  at  80  hz  provided  sufficient  reduction  in  the  high 
frequency  components  of  the  signal  that  the  feedback  gain  could  be 
increased  enough  to  significantly  affect  the  dynamics.  Unfortunately 
this  shifted  the  phase  of  the  feedback  signal  90°  at  the  fundamental 
bending  modes,  which  is  effectively  velocity  feedback  (i.e.  a  viscous 
damper).  The  compromise  was  to  implement  only  one  low-pass  filter  with 


B  -  1 


a  break  point  of  50  Hz  in  the  acceleration  feedback  path.  With  this 
configuration  the  feedback  was  still  effectively  an  acceleration  signal 
while  the  magnitude  of  the  high  frequency  noise  was  somewhat  reduced. 

3)  Some  high  frequency  oscillation  is  ocassionally  still  observed 
in  the  system,  especially  when  all  the  actuators  are  on  the  beam  and 
high  frequency  vibration  of  one  will  be  measured  by  another.  This  is 
not  a  consistent  problem  and  usually  doesn't  affect  the  overall 
experiment.  The  origin  seems  to  be  the  actuator  150  hz  bracket  dynamics 
still  feeding  back  through  the  compensator.  The  oscillation  is  not 
strongly  unstable,  and  usually  is  not  noticable.  Occasionally,  the 
response  gets  large  enough  that  it  is  audible  as  a  low  hum  coming  from 
the  actuator.  When  this  occurs  it  will  usually  stop  by  holding  the 
offending  actuator  on  or  around  the  LVDT  mounting  bracket.  The  signal 
can  also  become  noticeable  if  the  actuator  proof  masses  are  allowed  to 
rest  against  the  stops  when  not  commanding  the  system. 

A)  The  actuator  proof  masses  will  tend  to  drift  over  on  a  long 
term  basis  during  operation  due  to  a  very  low  frequency  drift  in  the 
power  amplifier.  The  power  amplifer  has  been  observed  to  oscillated 
with  a  period  of  approximately  AO  minutes  and  a  peak  voltage  of  0.15  V 
after  turn  on.  This  effect  damps  out  with  time,  but  the  actuators  have 
still  required  centering  adjustment  even  after  several  hours  of 
operation.  The  cause  of  this  oscillation  was  not  isolated  since  its  was 
of  such  a  low  frequency  relative  to  the  time  required  to  run  a  test  on 
the  beam.  Possible  causes  lie  in  :he  power  amplifier  itself  and 
possibly  in  the  power  supply.  This  problem  will  cause  the  proof  masses 
to  drift  into  the  stops  if  left  unattended  for  several  minutes.  Thus, 
proof  mass  centering  be  checked  prior  to  running  any  tests  on  the  beam. 


5)  A  configuration  for  the  compensator  was  found  that  provided 
acceptable  freqeuncy  response  performance  over  the  desired  bandwidth. 
However,  the  actuator  exhibited  a  drift  with  input  amplitude.  As  the 
amplitude  of  the  input  command  was  increased,  the  zero  reference  of  the 
power  amplifier  appeared  to  change.  This  caused  the  proof  mass  to  shift 
coward  one  end  of  the  travel  and  oscillate  about  that  point.  The  effect 
was  not  noticable  in  low  frequency  response  since  the  proof  mass  stroke 
was  saturated  before  the  amplitude  was  high  enough  to  cause  the  problem. 
However,  the  problem  was  noticable  at  the  first  torsion  frequency  of 
the  beam.  A  test  run  using  a  combined  fundamental  and  torsion  freqeuncy 
command  determine  that  the  maximum  output  force  that  could  be  applied  to 
any  two  modes  simultaneously  was  0.1  Ibf.  This  low  level  of  force 
output  was  not  acceptable  if  the  closed  loop  control  result  were 
expected  to  increase  the  beam  damping  by  more  than  a  few  percent. 

Various  approaches  were  tried  to  alleviate  this  problem  such  as  AC 
coupling  the  input  signal  to  ensure  it  has  no  bias;  AC  coupling  the 
com.mand  from  the  compensator  to  the  power  amplifier;  adding  another 
feedback  loop  around  the  power  amplifier,  which  was  found  to  have  other 
benefits,  but  had  no  affect  on  the  drift  problem;  and  substantially 
increasing  the  position  feedback  gain,  to  the  detriment  of  the  frequency 
response.  The  only  change  chat  had  any  affect  on  the  problem  was  the 
increase  of  the  feedback  gain.  However,  all  it  really  did  was  to  raise 
the  point  at  which  the  drift  became  noticable,  and  it  completely 
destroyed  the  flat  frequency  response  in  the  low  end  of  the  bandw’idth. 

Since  the  problem  appeared  to  be  associated  with  the  power 
amplifier  circut,  Mr.  Ken  Taylor  of  TRW,  who  did  the  original  design  of 
the  power  amplifier  circuit,  was  consulted  as  to  the  problem.  Mr. 
Taylor's  recommendation  was  to  check  the  very  high  frequency  output  of 
the  power  amplifier.  He  felt  the  amplifier  may  have  unstable 
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oscillations  in  the  kilohertz  to  megahertz  frequency  range  due  the  high 
impedance  of  the  input,  and  the  fact  that  the  amplifier  was  driving  a 
large  inductive  load  (the  motor)  through  very  long  cables  (8  to  10  ft 
for  this  experiment) .  He  suggested  adding  an  RC  network  across  the 
amplifier  output  if  the  high  freqeuncy  oscillations  were  observed. 

Since  the  Ono  Sokki  Analyzer  is  limited  to  40  KHz,  the  amplifier 
output  was  checked  on  an  oscilloscope.  A  ±15  V  oscillation  was  observed 
at  approximately  66  KHz.  Applying  the  fix  that  Mr.  Taylor  suggested  was 
successful  in  eliminating  it.  This  modification  did  not  change  the 
frequency  response  of  the  actuators,  and  did  relieve  the  input  amplitude 
drift  problem  to  some  extent.  The  drift  was  still  present  at  high  input 
amplitudes,  but  the  high  frequency  output  was  increased  to  almost 
2.5  Ibf  from  the  0.1  Ibf  limit  previously  measured.  While  this  is  still 
below  the  actuators  theoretical  capability,  it  is  a  much  more  acceptable 
level  for  this  experiment.  Later  testing  revealed  in  fact  that  this 
drift  was  not  an  actuator  problem,  but  was  due  to  a  very  slight  shift  in 
the  input  bias  from  the  Ono  Sokki.  The  actuators  were  found  to  respond 
to  as  little  as  2  mV  of  bias,  and  5  mV  would  run  the  proof  mass  to  the 
end  of  its  travel.  This  sensitivity  to  bias  and  very  low  frequency 
inputs  continues  to  be  a  problem  and  requires  very  high  quality  input 
signal  to  keep  the  actuators  centered. 


6)  The  phase  response  of  the  transfer  function  was  observed  to 
change  with  input  amplitude.  Some  variations  with  amplitude  were  also 
observed,  however  they  were  much  smaller  than  the  phase  shifts.  The 
initial  thought  would  be  that  as  the  input  amplitude  is  reduced,  the 
force  applied  by  the  motor  to  the  proof  mass  is  reduced.  The  expected 
consequence  would  be  that  friction  starts  to  become  a  significant 
portion  of  the  force  on  the  proof  mass,  and  accounts  for  the  phase 
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change.  Unfortunately,  the  shift  in  phase  is  opposite  to  the  expected 
effect  of  friction.  One  would  expect  friction  to  cause  a  lag  in  the 
output  as  it  holds  the  proof  mass  back  slightly  at  the  end  of  each 
stroke.  The  observed  behavior  however,  is  that  the  phase  increases  as 
the  input  magnitude  is  decreased.  In  an  attempt  to  determine  if  the 
phase  shift  was  related  to  friction  non-linearities,  the  time  response 
was  plotted  to  see  if  the  output  had  any  discontinuities.  Figures  B.l 
through  B.3  show  the  proof  mass  position  as  measured  by  the  LVDT 
compared  with  the  input  command.  The  data  indicates  that  there  is  no 
observable  sticking  of  the  proof  mass,  even  at  very  low  input  levels. 

This  problem  was  not  resolved,  but  the  phase  shift  from  zero 
degrees  in  the  final  configuration  is  within  ±10°  for  the  expected  input 
amplitude  range.  Isolation  of  this  problem  would  be  a  topic  to 
investigate  to  improve  the  accuracy  of  the  experiment. 

7)  hearing  and  shaft  friction  was  found  to  cause  a  problem  when 
the  actuators  were  mounted  on  the  beam.  The  proof  masses  float  to  a 
great  extent  as  the  beam  vibrates  but  still  add  approximately  3.5% 
damping  to  the  fundamental  bending  modes.  The  displacement  eventually 
damps  down  to  a  small  enough  level  that  the  acceleration  feedback 
command  cannot  overcome  the  friction  between  the  bearings  and  shaft.  At 
this  point  the  proof  masses  stop  floating  and  become  part  of  the 
structural  mass.  The  result  is  a  nonlinear  structure.  In  examining  the 
actuator  shafts,  scoring  by  the  bearings  was  becoming  significant.  The 
problem  became  acute  enough  in  actuator  #2  that  the  proof  mass  would 
occassionally  stick  while  being  commanded  and  often  required  manually 
moving  the  proof  mass  to  restore  free  motion.  At  this  time  the  problem 
still  exists.  New  actuator  shafts  with  a  much  harder  finish  are  under 
manufacture.  These  should  solve  the  scoring  problem,  and  may  help  the 


Response  -  High  Input  Level 


*■^51  w  <M  .w  m'  VA  w  ^ 


Fip.urp  B.2,  Proof  Mass  Position  Time  Response  -  Intermediate  Input  Level 


Respons 


low  amplitude  sticking.  Another  modification  that  will  further  reduce 
the  amplitude  at  which  the  proof  masses  stop  floating  is  to  increase  the 
acceleration  feedback  gain.  At  this  time  the  gain  is  as  high  as  is 
practical  for  this  configuration.  The  limiting  factor  in  the 
acceleration  feedback  has  been  high  frequency  noise.  The  feedback 
circuit  currently  has  one  low-pass  filter  on  the  accelerometer  signal. 
More  signal  conditioning  will  be  required  to  keep  high  frequency 
oscillations  out  of  the  actuator  when  the  gain  is  increased.  Increased 
acceleration  gain  will  also  affect  the  actuator  dynamics,  moving  the 
complex  pole  pair  to  even  lower  frequency.  This  may  actually  be 
beneficial  if  the  resonance  can  be  moved  below  the  range  of  any  low 
frequency  drift  that  might  be  in  the  command  signal.  If  this 
modification  is  made,  the  acceleration  feedback  amplifier  output  should 
be  checked  to  ensure  the  signal  is  not  close  to  saturating  either  the 
feedback  amplifier  or  the  subsequent  summing  junction. 


Appendix  C:  Actuator  Mass  Model 


To  develop  an  accurate  actuator  mass  model  the  configuration  was 
divided  into  six  components,  each  of  which  could  be  approximated  by 
simple  figures  (see  Figure  C.l).  These  components  include  all  of  the 
major  parts  of  the  actuator  system  with  significant  mass.  The  mass  of 
the  remaining  parts  not  represented  directly  by  one  of  the  components 
was  applied  to  the  nearest  representative  figure.  A  description  of  each 
component,  its  mass  and  dimensions  is  listed  in  Table  C.l. 


Table  C.l.  Actuator  Mass  Model  Components 


Figure 

Description 

Dimensions  (in) 

Length 

Width 

Height 

Diameter 

1 

Rectangular 

Plate 

142 . 34E-6 

1.1 

0.3 

1.7 

- 

2 

Circular 

Cylinder 

199 . 28E-6 

B 

- 

- 

0.25 

3 

Rectangular 

Plate 

483.93E-6 

1 

1.1 

0.3 

1  .  7 

- 

u 

Rectangular 

Plate 

199.28E-6 

1.1 

0.3 

2  .125 

• 

5 

Circular 

Cyl inder 

540.89E-6 

5.3 

J 

- 1 

0.75 

6 

Hoi  low 
Circular 
Cylinder 

2.9324E-3 

2.25 

Inside  Diameter  =  0.25 

Outside  Diameter  ^  2.125 

The  mass  moments  of  inertia  about  the  component's  center  of  gravity 
are  calculated  by  assuming  the  center  of  gravity  lies  at  the  center  of 
figure  then  orienting  parts  in  the  global  coordinate  frame  of  the 


structure  according  to  the  orientation  of  the  specific  actuator.  The 
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The  modal  test  of  the  final  experiment  configuration  was  conducted 
in  two  phases.  First,  the  fundmental  bending  modes  in  the  z-axis  and 
y-axis  were  measured  using  a  sine  dwell  technique.  Magnets  were  mounted 
on  either  side  of  the  base  plate  in  the  plane  of  the  bending  axis  to  be 
tested.  Magnetic  induction  coils  were  used  in  conjunction  with  these 
magnets  to  apply  a  sinusoidal  forcing  function  to  the  structure.  The 
frequency  and  damping  of  the  modes  were  determined  by  generating  a 
frequency  response  function  for  the  structure  around  the  natural 
frequency.  The  beam  was  forced  at  a  single  frequency  and  the  steady 
state  acceleration  of  the  base  plate  recorded.  Measuring  the  response 
at  intervals  throughout  a  frequency  band  around  the  being  mode  natural 
frequency  creates  the  response  function.  The  peak  response  occurs  at 
the  modal  frequency  and  the  damping  can  be  determined  using  the 
half -power  point  method  (equation  8).  To  determine  the  mode  shape,  the 
stiuctu’-c  was  forced  at  the  natural  frequency  and  the  acceleration  of 
eac-h  of  the  nodal  points  measured.  For  steady  state  motion,  the 
displacement  is  proportional  to  acceleration  and  is  therefore  a  direct 
irrasurement  of  the  mode  shape.  This  mode  shape  can  be  normalized  to  a 
uni;  vector  for  comparison  with  the  finite  element  data.  The  frequenev 
It  , npon.se  functions  are  shown  in  Figures  D.l  and  D.2,  and  the  measured 
modal  amplitudes  are  listed  in  Table  D.l. 

The  second  phase  used  the  GenRad  Modal  Analyzer  to  measure  the 
first  torsion  and  higher  bending  modes.  The  test  was  conducted  using 
thf-  impact  method.  The  beam  is  tapped  with  a  force  hammer  at  a 
selected  node  and  the  response  at  each  node  is  recorded.  The  measured 
time  response  is  transformed  to  the  frequency  domain  using  Fast  Fourier 
Transforms,  and  a  frequency  response  function  is  generated  for  each  node 


point.  The  frequency  and  damping  of  the  mode  are  found  from  the  curve 
fit  parameters  for  the  response  function,  and  the  mode  shapes  are 
determined  by  the  magnitude  and  phase  of  the  response  at  each  node 
point.  Figure  D.3  shows  the  node  point  geometry  used  internally  in  the 
GenRad.  Figure  D.A  is  the  first  torsion  mode  curve  fit  and  mod-' 
parameters.  Figures  D.5  through  D.8  show  the  higher  bending  mode  modal 
parameters  and  mode  shapes. 
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Table  D.l.  Z  and  Y  Axes  First  Bending  Modes 


Node 

First  Z  Bending 
Modal  Amplitudes 

First  Y  Bending 
Modal  Amplitudes 

Acceleration  (g) 

Acceleration  (g) 

1 

0.000000 

0.000000 

2 

0.000113 

0.000073 

3 

0.000252 

0.000674 

4 

0.000677 

0.001836 

0.001382 

0.003171 

6 

0.002461 

0.004519 

7 

0.003922 

0.006179 

8 

0.004738 

0.007805 

9 

0.007174 

0.010108 

10 

0.008762 

0.011995 

0.010791 


0.016117 
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3 
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4 

280.008 

0.000000 

6.008008 

5 

350.000 

8.000000 

0.000008 

6 

429.088 

0.000000 

0.000000 

7 

430.000 

0.008000 

0.000060 

8 

560.000 

0  000800 

0.008800 

<i 

630  000 

0.000000 

0  008888 

16 

8.000800 

0.000000 

0.006600 

14 

710.000 

0.008008 

0.000000 

Figure  D.3.  Modal  Test  Node  Point  Geometry 
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Figure  D.6.  Third  Z-Axis  Bending  Mode 
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Mod*  Sh«p«.  2ND  V-BENOINC 
Log  X  Coeff  Y  Coeff  Z  Coeff 

1  B  iddE-Bl  -b.?69E-«l  0  B80E-B1 

2  0  000E-81  -2.127E+00  0  000E-81 

3  0  000E-01  -3.825E+00  0  000E-01 

4  0  BeeE-ei  -s.sExE^Be  0  00eE-0t 

5  0  0e0E-ei  -6.684E+08  0  000E-ei 
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Figure  D.7,  Second  Y-Axis  Bending  Mode 
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Figure  D.8.  Third  Y-Axis  Bending  Mode 


Appendix  E:  Sensor  Specifications 


SPECIFICATIONS 

8628AS0 

8630A50 

8628AS 

*8630A5 

8692A50 

8692A5 

Acceleration  range 

g 

♦50 

♦5 

Acceleration  limit 

g 

♦80 

iS 

Sensitivity, ±5%  at  100  Hr,  3  g  nns 

mV/g 

100 

1000 

Threshold 

iiV  rms 

50 

100 

ug  rms 

500 

100 

Shock  (0.2nis  pulse  width),  max. 

g 

5000 

5000 

Amplitude  Linearity,  nom. 

% 

+  1 

*  ^ 

Time  Constant,  nom. 

s 

0.5 

0.5 

Resonant  Frequency,  mounted,  nom. 

kHz 

9 

9 

Frequency  Response, iS%  points 

Hz 

1...2k  1... 

Phase  Shift,  1...2k  Hz 

• 

<5 

<5 

Transverse  Sensitivity 

% 

<3 

<3 

Strain  Sensitivity  t  250u  in/in 

g/u  in/ in 

<.001 

A 

O 

o 

Temperature  Range,  operating 

•f 

+  32.  . 

.150 

Temperature  Range,  storage 

•f 

-10.  . 

.200 

Temperature  Coefficient  of  Sensitivity 

%/*F 

-.04 

-.04 

Supply  Current 

mA 

2. ..18  2... 

Source  Voltage 

VDC 

20.  .  . 

30  20.. 

Source  Impedance 

kfl 

100 

100 

Output  Bias,  nominal 

VDC 

11 

11 

Output  Impedance 

Q 

<100 

<100 

Output  Current 

mA 

2 

2 

Output  Voltage,  F.S.,  nominal 

V 

+  5 

♦5 

Ground  Isolation 

Material 

MR 

10 

10 

Housing  and  Base 

Aluminum  Alloy 

Seal 

Case  and  Connector 

Environmental 

Epoxy 

MODEL  8628  MODEL  8630 

MODEL 

8692 

weight  gr  6.7 
Mounting  Torque  in-lbf  8 


Figure  E,l.  Kistler  Accelerometer  Specifications 


Appendix  F:  Advanced  Beam  State  Space  Fortnulation 

Using  the  final  modal  test  data,  the  A  matrix  for  the  first  five 
modes  of  the  structure  is 


0 

0 

0 

0 

0 

-69.83 


0  -120.9 


0  -6247 


0  -12925 


0  -0.630  0 

0  0  -0.770 

0  0  0 

0  0  0 

•22545  0  0 


0  0  0 

1  0  0 

0  1  0 

0  0  1 

0  0  0 

0  0  0 

770  0  0 

0  -2.008  0 
0  0  -0,660 


0  -0.9610 


Applying  the  transfer  function  of  the  actuator,  the  ‘t>^B  matrix 
:akes  the  form 


<1>  B  = 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

2 . 747 

0 

-2.747 

-2.323 

0 

-2.323 

0 

0.484 

0.484 

0.496 

0.471 

0 

-0  215 

0 

0.215 

0.243 

0 

-0.243 

0 

Since  the  sensor  scale  factors  have  been  adjusted  to  vinitv  in  the 
input  to  the  PC-1000,  the  C  matrix  partitions  are  just  =  0,  a  4x5 
cero  matrix  and  C  has  the  form 

V 


c  <t>  = 

V 


0 

-4.943 

1.240 

0 

0.608 

4.905 

0 

1.240 

-0.573 

0 

0 

4.943 

1.240 

0 

-0.608 

-4.905 

0 

1.240 

0.573 

0 

These  matrices  can  be  partitioned  to  form  any  subset  of  modes  for 
inclusion  in  the  controller.  Examples  of  two  z-axis  controllers  are 
.sliown  in  Figvires  F.l  and  F.2.  These  configurations  include  the  first 
three  modes,  but  apply  control  gain  only  to  the  first  z  bending  mode. 
Estimator  and  controller  eigenvalues  differ  for  the  two  formulations  for 
cnirparison  of  feedback  gains. 
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Figure  F.l.  Z-Axis  Controller  Example  1 


::KIIV«.i£E  3F  C  •  KS  SrSTD!  ; 


l.uit;  « iMss 

I.U112  0.197333 

w.e«:  i.30oi2*-oi 

79  0415  1.20O12»-« 

lO  993*  V5370i=-C7 

10  9959  3.53Mle-02 


9.71896  6. 367999 
9.21890  0  302VH 
79.8105  0  102852 
^  |ll»  P,IW:S2 
11  2257  0.2(8(U 
11  7257  0.29|8» 


C.OOOO*«CO  C.00a0*»00  J.0900**0f  C.!50a?-^0!  0.C0C'JP*00  C.OOOOfOC 

:.ooco»*ca  '..oooo«*oc  c.ooeof+ca  o.ocooe«oc  o.‘/:50t*c;  c.ooooj'CC 

o.oooa»*oo  o.aoooe-ao  c.ooDte-*;::  i.oooc-e^o':  c.-je:;'?-:-.'  c.icace*;'; 

-C.£5E7(*‘'  ;.aCCO?»CC  C.??55?-'ar‘-c.i3;;e*o;  C.OoaaF*''C 

C.;003tfa0-J.;:j9f>j2  v.l'OSue-'S.'  v.ODCjp^iC 

0.0000»*CO  C.OOOO»fOC-C.t2<8*«W  O.OOOOt*aC  C.OOOOf‘OC-C.1897e»C; 

J 

0.0000r>00  0. 0000**00  0.0000**00  0.0000**00 
0.0000**00  0.0000**00  C.0000«*00  0.0000**00 

c.0000e*00  0.0000**00  C.0000»*00  0.0000**00 

0,0000**00  0.2747**0:  0.0000**0?-C.2747**0: 

-o.:::j*-3'.  o.ocoo*-*oo  c.;:;3**ci  ;. 0000**00 
v.<°:;»*00  0.4E:ie-OC  0.<96yt*0O  f..<7;ie*CC 


COCOe-OC- 

0000**70  0 

AAA*.g*.y  y 

OOCOs^O'-C 

CfO'je*?'*  C 

ocooe-cc  : 

00c’';t*07 

<?C®e*C'.  C 

i:;0*-'ic : 

A.'-.y.,*y  ' 

AAA/.g^.'*  A 

ocoo**C'  :• 

oooifo:  c 

OCOff-CC-'. 

4905**0’.  C 

oc.'0?*o;  c 

AAA.'.g*  *  *.« 

.oocoe-:  : 

. iTi-e-*'.  C 

.vv'.Ce*';' 

.::ne-':  : 

S  1 

v.'KKve-M)*' 

.292®e*0>C 

•  »  •  ^  -  V  V  .  ' 

."29e’rj'J 

O.i'T’e-OS  0.96O3e-0t-0.r9;6e-0T  j-Crtf-C'  c.;5:;t-c'-o.:;?5e-: 
■c, ;:76f-oi  o.4T::*-'.c-c.450i*-o5-o.49:;**oc  o.K75f-;:-t.. i:75t-c: 
C.90C2e-05-3.?60Ee-a5-0.i392e-C'7  0.4";e-C''’-C. 

:.::r4e-c:  c.4«:e-:c  c-.«otf-o5  ;.49:;**oc'  o.:r2*-;:-c.:2<:e-'. 


?.99E3i-o:  ;.99:3i-::  :.4;3Ec-::  :,£:sc:-;' 

-<.BS?::-'.9  9.?973:-r;;  7.;599;-;;  t.9“20-'.7  ;.:i'.0:-.‘.;  .".C:- 

;.;902;-'.o  ;.Bi45i-t;  9.B740;-o'.  -•..2:35;-ob  -3.”:ta-i:  i.c:'!?:- 

38450-01 -:.4l660-'.7  4.:B38:-9E  ?.e::9i-C:  :.::i3c-lt  -l.'oeO:- 

-i.::45a-;5-:,4C37d-.7;  9.40£5;-0E  ?.:470i-:<  9.85'.;:-i; 

5.c:;5:-o:  4.:r55-o9 -;.:24ie*o:  -i.iesEa-c:  -:.i536s-0E 


4. :585d-14  <.3C15d-C5  4.:5e3d-l*  -;.::;5:-75 
-:.4594d-C5  3.3?B0t-14  :.t5''4e-''5  Z.Vl'Ji-'.'- 

:.4SC5a-06  :.4505d-04  :.i505d-04  :.43C£a-i4 

5. :r:e-;; -:.oo49c-o:  i.zzzu-::  ;.)6*9t-o: 

;.9r:*:-o’  4.44;:d-:;  -i.crrt-o:  4.4«:;:-;; 
•5.7913d-9:  -«.79:3e-o:  -5.79;3:-;:  -5.T9;3c-t: 


Figure  F.2,  Z-Axis  Cont  rol 'ei  P.> 
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